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ABSTRACT 


Theoretical  studies  concerning  the  mechanics  of  penetration  and 
perforation  of  a  snow  covered  Arctic  sea  ice  sheet  subjected  to 
projectile  impact  were  performed. 

Penetration  problems  were  treated  by  a  deep  penetration  theory 
based  on  dynamic  spherical  cavity  expansion  analysis.  In  particular, 
finite  compactibility  and  permanent  deformation  of  both  the  snow  and 
sea  ice  target  materials  were  taken  into  account  by  assuming  a  locking 
approximation  for  behavior  under  hydrostatic  stress,  and  response  as  an 
elastic-plastic,  linear  strain  hardening  solid  under  shear  stress, 

Perforation  problems  were  treated  with  the  aid  of  a  two-dimensional, 
large  deformation,  dynamic,  elastic-plastic  computer  code  (CANDIA  CODE) 
which  was  developed  in  a  previous  investigation.  Specifically,  axi- 
symmetric,  dynamic  stress  distributions  were  studied  under  conditions  of 
impact  at  normal  incidence  for  a  cylindrical  blunt  end  projectile  and  a 
sea  ice  target  slab.  A  capability  for  considering  perforation  problems 
where  fracture  in  the  sea  ice  target  material  occurs  was  developed.  In 
this  connection,  the  redistributions  of  both  tensile  and  shear  stresses 
that  accompany  propagating  fracture  surfaces  were  acknowledged.  The 
program  was  also  modified  to  provide  projectile  deceleration  loads 
during  the  perforation  process,  and  subroutines  were  created  to  treat 
multilayer  ice  slabs  and  thereby  accommodate  effects  contributed  by 
overlying  snow  cover  and  the  underlying  mushy  sea  ice  skeleton  layer. 
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PREFACE 


The  investigation  reported  here  was  sponsored  by  the  U.S.  Naval 
Ordnance  Laboratory  (NQL,  White  Oak),  Silver  Spring,  Maryland,  under 
ONR  Contracts  Nonr-2332(00)  and  N00014-68-A-0243  (SRI  Project  No.  7000-452). 
The  work  was  performed  between  1  July  1969  -  31  December  1969  and 

15  April  1970  -  31  June  1970  and  is  designated,  unofficially,  Phase  IV 

( APAUXj )  (.!■  '  /  /<>• 

Study.  Results  of  the  'Phase  I,  Phase  II,  and  Phase  III  studies  were 
published  in  separate  Stanford  Research  Institute  reports  under  a 
similar  title  and  dated  November  1965,  May  1967,  and  September  1969. 

The  contract  was  monitored  for  the  Naval  Ordnance  Laboratory  by 
Mr.  M.  M.  Kleinerman.  Project  leader  for  Stanford  Research  Institute 
was  Dr.  Bernard  Ross,  and  principal  investigators  were  Dr.  Sathya  Hanagud, 
Dr.  Bernard  Ross,  and  Mr.  Gursharan  Sidhu. 
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12  Eq.  Ill -34 
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C  -  constant  of  integration 
o 

C  ,  C  -  constants  related  to  dynamic  pressure,  see 
1  2  Eq.  II I -77 

D  -  diameter  of  projectile 
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E  -  tangent  modulus  for  linear  strain-hardening 
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p(t)  -  dynamic  pressure  applied  to  spherical  cavity 
surface  as  a  function  of  time 
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point,  q  =  R 
°s 
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R  -  depth  of  snow-ice  interface  at  transition 
point  of  projectile  penetration  process, 
see  Fig.  2 
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r  -  radial  coordinate 

r'  -  dummy  variable  of  integration 
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V  -  velocity  of  projectile  after  completion  of 
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I  INTRODUCTION 


This  report  presents  the  final  results  of  a  comprehensive,  long 
term  investigation  concerning  the  penetration  and  perforation*  behavior 
of  Arctic  sea  ice  subjected  to  impact  by  inert  projectiles.  The  impor¬ 
tance  of  these  engineering  problems  to  the  Arctic  Antisubmarine  Warfare 
Program  has  been  described  in  three  reports  under  the  present  title 
issued  during  November  1965,  **  May  1967,  and  September  1969.  The 
last  report  contained  a  complete  summary  of  work  performed  and  progress 
achieved  up  to  that  date;  therefore,  the  topic  is  not  considered  anew. 

However,  it  is  noted  briefly  that  two  important  analytical 
procedures  for  treating  projectile-sea  ice  interaction  problems  were 
developed  after  completion  of  the  second  study.  As  described  in  the 
last  report,  these  efforts  included  the  realization  of  a  large  defor¬ 
mation,  elastic-plastic,  artificial  viscosity  type  computer  program 

(CANDIA  CODE)  which  could  be  used  to  solve  a  wide  variety  of  axisymmetric , 

4  5  6  7  8 

two-dimensional,  dynamic,  and/or  impact  problems,  ’  ’  ’  ’  and  a  large 

deformation  deep  penetration  theory  for  analyzing  projectile  motion  in 

. .  9,10,11,12 

compactible  media. 

Even  though  it  was  possible  to  derive  a  substantial  body  of 
information  by  applying  these  procedures  to  Arctic  sea  ice  perforation 


*  For  the  purpose  of  this  study,  penetration  can  be  defined  as  'he 
entrance  of  a  projectile  into  a  target  or  ice  cover  without  com¬ 
pleting  its  passage  through  the  body,  whereas  perforation  implies  the 
complete  piercing  of  the  target  slab  or  ice  cover  by  the  projectile. 

**  Superscript  numbers  refer  to  references  which  are  collected  at  the 
end  of  this  report . 
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and  penetration  problems,  it  was  recognized  that  further  improvements 
were  still  required  for  practical  application  of  the  theoretical  work. 

For  example,  perforation  of  an  Arctic  sea  ice  cover  by  an  impacting 
projectile  creates  fracture  surfaces  in  the  target  material,  and  it  is 
the  resultant  geometry  of  this  process  that  describes  both  the  mode  of 
perforation  and  the  associated  value  of  critical  impact  velocity.* 
Moreover,  important  redistributions  of  stresses  occur  when  fracture 
surfaces  develop,  and  these  changes  affect  both  the  failure  phenomenon 
itself  and  the  magnitudes  of  stress  experienced  by  the  Impacting  pro¬ 
jectile.  In  this  connection,  the  previous  version  of  the  code  was 
capable  of  providing  complete  descriptions  of  both  stress  and  defor¬ 
mation  fields  in  the  impacting  bodies.  It  remained  to  postulate  suitable 
failure  criteria  and  expected  crack  propagation  paths  to  define  the 
global  fracture  pattern  and  the  mechanism  of  projectile  transit  through 
the  sea  ice  cover. 

In  consequence  of  the  present  work,  it  is  possible  now  to  track  the 
actual  failure  configuration  as  a  function  of  time  after  initial  impact, 
and  to  acknowledge  that  the  complement  of  both  shear  and  cleavage  frac¬ 
ture  surface  propagation  paths  is  being  influenced  continuously  by  the 
developing  fracture  process. 

On  the  other  hand,  the  previous  deep  penetration  theory  was  limited 
in  application  to  homogeneous,  isotropic  media  and  thus  was  not  suitable 
for  cases  involving  a  snow  covered  Arctic  ice  sheet.  To  correct  this 
shortcoming,  present  research  efforts  were  directed  toward  realization  of 
a  multilayer  deep  penetration  theory  to  treat  the  practical  projectile 

*  Critical  velocity  is  defined  as  the  minimum  or  threshold  impact 
velocity  for  complete  penetration,  i.e.,  perforation. 
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penetration  problem  characterized  by  a  two-layer,  snow-ice  target 
material.  This  work  can  be  employed  to  obtain  penetration  performance 
curves  for  projectiles  deployed  in  a  typical  Arctic  environment. 

In  sum,  these  connected  analyses  (namely,  direct  treatment  of  the 
dynamic  fracture  problem  and  synthesis  of  a  multilayer,  large  deformation, 
deep  penetration  theory)  constitute  the  basis  of  this  investigation. 
However,  in  addition  to  these  studies  a  subsidiary  effort  concerning  the 
determination  of  maximum  deceleration  loads  experienced  by  a  projectile 
during  the  sea  ice  perforation  process  was  carried  out.  Unfortunately, 
numerical  calculations  for  both  of  these  problems  could  not  be  carried 
out  due  to  a  lack  of  funds  for  computer  runs.  Therefore,  the  significance 
of  present  developments  as  applied  to  the  design  and  development  of  naval 
weapons  for  Arctic  application  could  not  be  examined. 


II  MECHANICS  OF  FRACTURE 


The  work  accomplished  was  concentrated  on  modifying  and  restructuring 
the  CANDIA  CODE  to  enable  theoretical  analysis  of  problems  where  fracture 
of  the  ice  slab  occurs.  This  phenomenon  engenders  redistributions  of 
stresses  in  the  target  material,  and  accompanying  effects  in  the  magni¬ 
tude  and  distributions  of  loads  sustained  by  the  impacting  projectile. 
Thus,  from  an  overall  viewpoint,  these  efforts  were  directed  toward  the 
transformation  of  a  computer  code  employed  previously  as  a  research  tool 
to  a  working  code  which  can  be  utilized  to  perform  engineering  design 
and  development  calculations. 

To  sustain  these  particular  goals  and  to  treat  the  fracture  problem 

in  sea  ice,  per  se,  it  was  necessary  to  recast  the  code  so  that  problems 

characterized  by  greater  numbers  of  variables  and  finer  finite-difference 

calculation  mesh  networks  could  be  accommodated.  For  example,  in  initial 

format  the  code  considered  sea  ice  to  be  homogeneous  and  isotropic, 

whereas,  in  fact,  this  material  is  strongly  anisotropic  and  inhomogeneous. 

Consequently,  introduction  of  this  more  advanced  description  of  material 

composition  would  furnish  a  demand  for  increase  in  the  number  of  input 

variables  that  are  needed  for  numerical  calculation.  Moreover,  greater 

machine  storage  and  capability  were  also  required  because  a  developing 

fracture  surface  submits  many  more  mesh  points  of  importance  into  the 

calculation  process.  To  overcome  these  difficulties,  the  entire  CANDIA 

3 

CODE,  as  listed  and  flow-charted  in  the  last  comprehensive  report,  was 
restructured  to  eliminate  its  dependence  on  high  speed  core  memory 
storage.  As  a  result,  the  code  in  its  present  format  can  be  employed  to 
take  advantage  of  random  access  disc  storage. 


Preceding  page  blank 


5 


At  this  time,  the  entire  code  exists  on  disc  in  precompi led  form 
which  is  compatible  with  all  electronic  data  processing  equipment  based 
on  the  FORTRAN  language.  As  such,  the  code  gains  advantage  in  capabil¬ 
ity  and  versatility  by  employing  fixed  or  consistent  amounts  of  high 
speed  core  storage  with  a  random  access  interface. 

Briefly,  advantages  of  the  new  disc  format  are: 

•  Capacity  increased;  that  is,  problems  characterized 
by  greater  numbers  of  intrinsic  variables  and  re¬ 
quiring  finer  mesh  works  for  calculation  purposes 
can  be  treated.  The  fracture  problem  in  sea  ice 
falls  within  this  definition. 

•  Convenience  Increased;  that  is,  ease  of  working 
with  the  code  has  been  considered.  For  instance, 
the  code  is  transcribed  presently  on  one  disc  pack 
in  lieu  of  multiple  racks  of  perforated  cards. 

•  Versatility  increased;  that  is,  interrupted 
calculations  can  be  performed  so  that  when  a 
computation  procedure  is  stopped  intermittently 
it  is  no  longer  necessary  to  resume  calculations 
from  the  initial,  or  zero,  input  state.  Instead, 
numerical  treatment  can  be  restarted  from  the  pre¬ 
vious  termination  point.  In  this  manner,  and  with 
particular  relevance  to  the  fracture  problem,  an 
analysis  can  be  intercepted  to  focus  attention  on 
the  development  of  a  propagating  fracture  surface 
and  the  attendant  stresses  in  plastic  zones 
surrounding  a  crack  tip. 

•  Information  display  possibilities  increased;  that 
is,  direct  computer  readout  is  feasible  to  furnish 
continuing  computer  code  output  by  means  of  cathode 
ray  terminals.  Thus,  it  becomes  practical  to  take 
movies  of  progressive  display  patterns  and  thereby 
indicate  visually  the  continuing  propagation  of 
fracture  surfaces  and  ancillary  stress  redistri¬ 
butions  in  both  target  material  and  impacting 
projectile. 

With  restructuring  and  modification  of  the  computer  code  effected, 
it  was  possible  to  attack  the  fracture  problem  directly.  However,  in 
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this  case  it  was  acknowledged  first  that  the  treatment  of  stress  redis¬ 
tributions  due  to  developing  fracture  surfaces  would  be  closely  coupled 
to  the  states  of  stress  existing  on  free  surface  boundaries  such  that 
the  proper  numerical  treatment  of  boundary  conditions  themselves  would 
be  an  utmost  necessity.  In  this  connection,  it  was  pointed  out  in  Ref.  3 
that  many  large  deformation  elastic-plastic  computer  codes  currently  in 
existence  were  founded  on  incorrect  boundary  condition  formulations. 

Thus,  an  appendix  was  included  in  the  report  containing  a  lengthy  deri¬ 
vation  for  the  true  treatment  of  these  boundary  conditions.  Because 
this  development  was  somewhat  unwieldy  in  the  form  presented,  consider¬ 
able  effort  was  expended  during  the  current  period  in  rewriting  a  more 
concise  and  simplified  algorithm  for  this  aspect  of  the  problem. 

Once  this  task  was  completed,  it  was  possible  to  create  an 
additional  algorithm  to  treat  the  fracture  problem  itself.  In  this 
case,  the  algorithm  was  founded  on  the  development  of  a  point  condition 
code  which  onabled  integer  point  condition  values  to  be  associated  with 
different  characteristic  free  surface  points  on  nodal  intersections  in 
the  finite  difference  network.  For  example,  individual  point  values 
were  ascribed  with  concomitant  subroutines  to  represent  lateral  and 
transverse  free  surface  points,  interior  fracture  surface  points,  and 
corner  points  on  both  target  and  projectile.  As  a  result,  point  con¬ 
dition  changes  could  be  related  to  responsive  subroutines  in  correspon¬ 
dence  with  physical  transformations  of  both  impacting  bodies.  Conse¬ 
quently,  it  is  possible  now  to  track  in  the  computer  both  current 
boundary  states  and  conditions,  as  well  as  emerging  and  propagating 
fracture  surfaces. 

Use  was  then  made  of  the  modified  boundary  condition  and  fracture 
algorithms  to  treat  the  shear  fracture  problem  encountered  when  a  blunt 
or  hemispherical  nose  projectile  experiences  initial  penetration  in  a 
target  material.  In  this  case,  it  was  known  from  previous  experimental, 
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field  test,  and  theoretical  studies  that,  due  to  the  anisotropic  nature 
of  a  sea  ice  target  material,  initial  fracture  under  these  impact  con¬ 
ditions  would  occur  in  the  shear  mode  along  c  -normal  stress  plsnes  in 

zz 

the  vertical  direction.* 


With  the  shear  fracture  phase  of  the  overall  perforation  process 
accounted  for,  work  progressed  toward  realization  of  a  fracture  algorithm 
to  describe  failure  surface  propagation  under  tensile  principal  stresses. 
The  Important  difficulty  overcome  in  treating  this  problem  was  that  frac¬ 
ture  does  not  occur  necessarily  along  a  mesh  direction.**  Therefore, 
there  remained  the  burden  of  outlining  and  defining  not  only  the  stress 
redistributions  at  incipient  fracture  but  the  directions  of  the  fracture 
planes  themselves. 


In  this  context,  a  method  of  approach  was  devised  and  implemented  in 


the  present  version  of  the  code. 


Specifically, 


all  components  of  the  a 


ij 


*  Fracture  configurations  in  the  shear  mode  are  initiated  upon  projectile 
impact  at  normal  incidence  because,  for  a  blunt  penetrator,  regions  of 
intense  shear  stress  are  generated  around  the  projectile  periphery. 
Moreover,  laboratory  experiments  and  field  test  results,  as  well  as 
theoretical  calculations  (CANDIA  CODE),  indicate  that  continued  shear 
penetration  and  concomitant  fracture  patterns  are  engendered  in  the 
sea  ice  slab  until  penetration  depths  of  the  order  of  70-80%  sea  ice 
sheet  thickness  are  achieved.  At  this  point,  it  is  recognized  that 
shear  stress  magnitudes  fall  off  and  that  further  fracture  and  ultimate 
perforation  are  the  result  of  sea  ice  failure  under  tensile  stress.  In 
sum,  complete  penetration  is  obtained  through  a  composite  process  char¬ 
acterized  by  initial  failure  due  to  shear  and  terminal  failure  as  a 
result  of  tensile  or  cleavage  fracture.  Consequently,  ultimate  perfo¬ 
ration  results  in  the  ejection  of  a  cylindrical-conical  shear  plug  of 
target  material  by  the  impacting  body. 


**  In  this  connection,  the  orientations  of  component  segments  of  a 
developing  tensile  fracture  surface  are  related  directly  to  the 
magnitudes  and  directions  of  local  principal  stresses,  which  them¬ 
selves  are  dependent  variables  of  the  problem  in  question.  Thus, 
developing  cleavage  fracture  surfaces  do  not  correspond  necessarily 
in  spatial  alignment  and  configuration  to  the  postulated  geometry  of 
the  mesh  network  employed  for  finite-difference  calculations. 
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stress  tensor  are  determined  at  the  center  point  of  a  given  area  segment 
in  a  mesh  network.  Knowing  these  quantities  enables  the  three  principal 
stresses  to  be  determined,  whereupon  a  test  is  performed  to  isolate  any 
possibility  that  the  shear  stress  component  of  the  total  stress  tensor 
violates  the  allowable  shear  failure  stress  of  the  target  material. 
Assuming  this  is  not  the  case  and  that  shear  fracture  does  not  develop, 
then  the  calculation  sequence  is  pursued  further  by  comparing  the  maximum 
principal  stress  value  with  the  allowable  failure  stress  in  uniaxial 
tension.  If  the  former  quantity  exceeds  the  latter,  it  can  be  assumed 
(i.e.,  at  least  for  an  elastic  brittle  material  such  as  freshwater  ice) 
that  fracture  occurs  over  the  mesh  section  of  interest.  Then,  spatial 
orientation  of  the  fracture  surface  segment  is  obtained  by  considering 
the  cleavage  plane  to  be  orthogonal  to  the  direction  of  the  maximum 
principal  stress  component. 

Thus,  knowledge  of  the  magnitudes  and  directions  of  principal  stress 
quantities  at  a  central  point  in  the  mesh  section  enables  the  prediction 
of  fracture  likelihood  and  resultant  local  failure  surface  geometry  to  be 
made.  Once  fracture  occurs  over  a  mesh  section,  it  is  necessary  to  pos¬ 
tulate  free  surface  boundary  conditions  which  provide  a  zero  traction 
state  over  the  developing  area  of  material  separation.  Then,  the  physi¬ 
cal  situation  at  the  failed  mesh  section  must  be  rechecked  at  each  and 
every  time  step  cf  succeeding  numerical  calculation  to  determine  whether 
changes  in  external  loading  and/or  internal  stress  state  have  resulted  in 
subjecting  the  mutual  fracture  surfaces  to  closure  and  compressive 
pressure.  If  this  turns  out  to  be  the  case,  free  surface  boundary  con¬ 
ditions  are  stricken  and  appropriate  interface  conditions  substituted  in 
their  place.  Finally,  at  the  conclusion  of  each  calculation  time  step, 
a  locus  or  contour  is  passed  through  all  of  the  individual  fracture  sur¬ 
face  segments  to  indicate  instantaneous  states  of  the  composite  failure 
pattern  developed  in  the  target  material  subjected  to  projectile  impact. 


9 


Ill 


MECHANICS  OF  PENETRATION 


In  addition  to  continuing  research  in  the  fracture  problem,  efforts 
were  devoted  to  realizing  a  theory  for  the  deep  penetration  of  projec¬ 
tiles  in  multilayered  media  such  as  a  snow  covered  Arctic  sea  ice 
13, 14,15 

sheet.  Solutions  of  the  dee  penetration  problem  for  a  single 

layer,  homogeneous,  isotropic  sea  ice  cover  have  been  obtained  already, 

and  their  formulation  and  development  were  presented  in  a  recently  issued 

3 

comprehensive  report.  Fortunately,  most  of  this  previous  work  was 
applicable  directly  to  the  multilayered  penetration  problem;  thus,  the 
compounding  of  a  tractable  theory  applicable  to  projectile-sea  ice  inter¬ 
action  problems  for  a  more  complicated  target  material  did  not  offer 
substantial  difficulty. 

In  this  case,  continued  use  was  made  of  the  two  most  important 
conceptual  fundamentals  upon  which  the  previous  theory  was  based.  These 
are:  first,  that  transformation  of  the  deep  penetration  problem  itself 
can  be  effected  successfully  to  yield  a  problem  concerned  with  dynamic 
stress  distributions  in  an  infinite  solid  that  contains  a  spherical 
cavity  subjected  to  step  input  in  pressure*,  and,  second,  that  the  dynamic 
constitutive  behavior  of  compactible  materials  such  as  Arctic  sea  ice 
and  snow  can  be  represented  satisfactorily  by  assuming  an  idealized 
locking  approximation  for  Rankine-Hugoniot  material  behavior  under 
hydrostatic  compressive  stress. 

*  This  metamorphosis  is  founded  specifically  on  the  assumption  that  dy¬ 
namic  pressure  experienced  by  a  penetrating  projectile  at  its  stagnation 
point  is  comparable  by  reason  of  symmetry  to  the  similar  pressure  quan¬ 
tity  which  exists  at  the  surface  of  a  dynamic  expanding  spherical  cavity. 
A  corollary  assumption  enables  an  expression  for  the  complete  force 
resisting  continued  penetration  to  be  established  by  postulating  that 
the  spatial  variation  of  pressure  over  the  frontal  portion  of  the 
projectile  can  be  represented  by  a  simple  cosine  relationship.1" 
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Preceding  page  blank 


Both  of  these  arguments  were  invoked  to  resolve  the  multilayer 
penetration  problem.  In  particular,  a  systematic  approach  was  devised  to 
recast  the  two-layer  target  configuration  into  the  form  of  an  idealized 
infinite  body  that  included  a  spherical  cavity  surrounded  by  two  concen¬ 
tric  spherical  regions  of  material  having  the  relevant  physical  and 
mechanical  properties  of  Arctic  snow  and  sea  ice.  Then,  the  dynamic 
spherical  expansion  problem  was  worked  through  again  for  the  composite 
material  configuration.  This  task  demanded  relatively  more  complicated 
algebraic  formulations  but  did  not  require  an  increase  in  mathematical 
generality  or  Important  theoretical  innovation.  The  problem  was  made 
somewhat  easier  by  the  fact  that  reflected  stress  waves  in  a  locked 
material  travel  at  infinite  speeds,  so  that  introduction  of  limiting 
conditions  into  the  mathematical  process  resulted  in  well-defined 
simplifications . 

The  problem  considered  is  depicted  schematically  in  Figs.  1  and  2. 

Attention  is  directed  first  to  the  ideal  locked-rigid  plastic  material 

in  Region  1  (Figs.  3  and  4)  which  is  behind  the  advancing  shock  front 
s 

at  r  =  b(t) .  The  locking  condition  on  strain  is 


e  +  2eQ  =  «,  (III-l) 

r  6  *s 

where  the  ideal  locking  strain  for  snow  is  a  time-invariant  material 
property.  Differentiating  with  respect  to  time  results  in 

C  +  2*e  =0  (III-2) 

r  e 

These  strain  rate  quantities  can  be  expressed  in  terms  of  radial  outward 
particle  velocity,  v(r,t),  by  the  kinematic  relations 


c 


r 


(Sv 

3r 


(III —3) 
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FIGURE 


SCHEMATIC  REPRESENTATION  OF  PROJECTILE 
PENETRATION  PROCESS  IN  SNOW  COVERED 
ICE  SHEET  INITIAL  ENTRY  AND  FIRST 
PHASE  MOTION 


FIGURE  2  SCHEMATIC  REPRESENTATION  OF  PROJECTILE 
PENETRATION  PROCESS  IN  SNOW  COVERED 
ICE  SHEET  TRANSITION  AND  SECOND  PHASE 
MOTION 
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SNOW 


DILATATION  SHEAR 

x  denotes  condition  of  MATERIAL  IN  REGION  If 

SNOW 


FIGURE  4  IDEALIZED  STRESS- STRAIN  CURVES  FOR 

A  RIGID  PLASTIC,  IDEAL  LOCKING  MATERIAL 
WITH  LINEAR  STRAIN  HARDENING  AFTER  VIELD 
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Substitution  and  integration  yield 


f  (t) 

v  =  —  (III  -4) 

s  r 


where  subscript  s  denotes  quantities  in  the  locked-plastic  snow 

region.  The  stress-strain  relationship  in  shear  for  a  rigid-plastic 

3 

material  with  linear  strain  hardening  is 

a.  -  a  =  Y  +  ~  £t  (€.  -  c  )  (III-5) 

8  r  s  3  *s  8  r 

Under  conditions  of  spherical  symmetry,  the  equation  of  motion  becomes 


'ov 

s 

l  z —  +  v 
s  \  ot  s 


( 1 1 1  -6 ) 


Finally,  the  large  deformation  circumferential  strain  is  given  by 


(III ~ 7 ) 


Combining  these  relationships  into  Eq.  Ill -6  and  integrating  provides 
the  following  expression. 


f2 

~7  +  g  (t) 
r  s 


(III -8) 


Conservation  of  mass  across  the  shock  front  at  the  stress  free-locked 
plastic  interface  yields 

pn  i)  =  p,  (i>  -  v  )  r  =  b(t)  (III-9) 

°s  *s  s 
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or  upon  reduction 


f  (t)  ■  or  b26 

S  S 


(III-10) 


where 


■  1  -  "o 


(Ill -1)) 


Similarly,  conservation  of  momentum  yields 


O  =  -Of  (b-v)v 

rS  *8  S  S 


r  =  b(t) 


(III-12) 


or  upon  reduction 


r=b 


PD  a 

°8  S 


b2 


(III-13) 


The  boundary  condition  on  stress  at  the  cavity  surface  is 


a_  =  -  p(t)  @  r  =  a(t) 


(III-14) 


Then,  from  Eq.  III-8  evaluated  at  r  =  a(t) 


,  2c  |  \ 

g  (t)  =  -  p(t)  -  2Y  In  a  -  2  /  ~  Et  1 2  In  —  +  —”)dr' 

s  /  r'  ls  V  r  3 

'f  N 

s 


a 

l* 


+  0 1 


s  1 


s  a  2  a 


(III -15) 
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so  that  the  equation  for  radial  stress  becomes 


+  "1.*.  (b2£  +  2bS>)  (a  '  r) 


1_ 

2 


n 


a  b4  b2 


f-L  _  -1\ 
Va4  r4  J 


(111-16) 


Using  Eq.  III-13,  an  expression  is  obtained  for  the  pressure  variation 
with  time  at  the  cavity  surface 


P(t) 


2Y 

s 


In 


+  2bb2 


Equation  III-17  can  be  reduced  further  by  application  of  the  boundary 
condition  on  particle  velocity  at  the  cavity  surface 

v  =  a  @  r  =  a(t)  (II I -18) 

s 

Using  Eqs.  Ill -4,  III -10  results  in 

a  b2  b  =  a2  a  (III -19) 

s 
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1  4  *a 

2  PV  8 


(IIX-21) 


Surfaces  of  discontinuity  in  displacement  and  particle  velocity  cannot 

exist  if  fracture  at  the  shock  front  is  omitted  from  consideration. 

Thus,  v  (t)  is  zero  on  the  shock  front  since  the  stress  free  material 
s 

ahead  of  the  locked  plastic  zone  is  quiescent.  As  a  consequence,  no 
material  moves  through  the  shock  front  and  the  following  compressibility 
relationship  can  be  applied  at  successive  times  over  the  volume  of 
locked  plastic  material  bounded  by  the  cavity  surface  and  an  arbitrary 
spherical  surface  at  the  present  shock  front  location 


Upon  reduction 


(III -22) 


(III-23) 
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16, 17, 18 

But,  it  has  teen  shown  by  Goodier  that,  for  the  deep  penetration 

problem,  the  following  assumption  can  be  made 


a3  »  a3 


In  addition,  f>,  >  P  ,  so  that 

jC  O 


b3  s  —  a3 
o 

s 


(III -24) 


Equation  III-21  becomes 


P(t)  = 


-  -  Y  In  cr  +2 
3  s  s 


+  P,s  («S  ♦  2iS)  (i  -  *») 


of'3)  -  o0  iV 

2*8  \  s  /  °s  s 


(III-25) 


The  remaining  integral  must  be  evaluated.  Application  of  an  appropriate 
compressibility  relationship  at  an  arbitrary  location,  r,  within  the 
locked  plastic  region  yields 


(r3  -  a3)pjts  =  (r3  -  a3)P0£ 


(III -26) 
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Using  the  deep  penetration  assumption  (i.e.,  a3  »  a3),  Eq.  III-26 
becomes* 


rc  i1/3 

°s  rVa3 

Oi  r3/a3  -  1 


(HI-27) 


Letting 


x  =  5s  =  r3/a3 


and  employing  Eq.  Ill -27  enables  the  integral  expression  in  Eq.  Ill -25 
to  be  written  as  follows 


(footnote  continued  on  next  page) 
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Evaluation  of  the  first  two  integrals  in  the  second  member  of  Eq.  Ill -28 

3 

is  straightforward.  The  remaining  integral  is  attacked  by  letting 
y  =  1/x  and  changing  the  lower  limit  of  integration  to  1 -e  where 
e  «  1.  There  results 

b3  /a3  o‘: .  bJ 


/ 


x  dx 
In  — -  — 
x-1  x 


/ 


In  { 1-y) 


dy 


l-€ 


(III -29) 


Since  b  >  a  always  and  1-s  <  1,  the  range  of  values  over  which  integration 
takes  place  (i.e.,  l~e  >  y  >  b3  a*3)  ensures  that  y  <  1  so  that  the  inte¬ 
grand  can  be  expanded  in  series  form.  Step  by  step  integration  yields 


b3/a3 


/ 


In 


dx 
x-1  x 


£_  £. 

+  +  —  +  ... 

4  9 


i  -e 


3  W3 
a  b 


(II I -30) 


(footnote  continued) 


but  p»  >  P„  and  r3  »  a3 
°s  °  o 


r3/rj  -  a3/r'0 


=  1 


'»„s/Pls)r 

5/3  3  /  5 

r  /rQ  -  a  /rQ 


=  r 


(p„  /p„  )r3r3 


°s 

~r3  -  a3 


-  r ' 


(0  /o»  )r*'/a1 

°s  <  s _ 

r  /'a'  -  1 
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Allowing  e  “►0  and  introducing  Eq.  Ill -24  results  in 


x  dx 

m  —  — 

x-1  x 


o  1  r 
-  >  t  a 
n  s 


(III-31) 


Finally,  Eq.  Ill -25  becomes 


2  4  rf*  x  „ 

-  ~  Y  In  a  +  -  Et  — —  )  ,  Ta 

3  s  s  9  ts  6  n  s 

n=l 

+  »1B  («*  ♦  2**)(l  -  <3) 

i1  ■  *,'*)  *  °o/“s1/3  <in-32) 


The  expression  for  pressure  can  be  broken  down  into  static  and  dynamic 


parts  so  that 


P(t)  =  Pg  +  Pls(\a*  +  b2®2) 


(III -33) 


where 


2  o  2  4  <C~''  l  n 

=  —  e*  -  -  Y  In  a  -  -  Et  >  —a 

27  ^s  3  s  s  9  ls  xr 


B  =  l  -  al/3 
l  s 


§  -  (l  +  a  )oiL/3  +  J  °P/: 

2  s  s  2  s 


-  -  (XL/3  -  -  a4/3 
2  s  2  s 


(III-34) 
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The  relationship  for  pressure  variation  at  the  cavity  surface  given  by 

3,16 

Eq.  III-33  is  used  in  the  Goodier  deep  penetration  theory  to  obtain 
an  expression  for  resisting  forces  acting  on  the  projectile.  Then,  an 
equation  of  motion  for  projectile  transit  in  the  target  material  can  be 
written 


Mq  =  "  PS  +  JP*s 


(®.  1 5 + b/)! 


rny 


(III-35) 


After  integration,  Eq.  III-35  becomes 


In  (PS  *  I  WO  =  '  3  B3°Js  5" 


7  +  B  -  0, 

A  i  3  *s 


+  C  (111-36) 
o 
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Exact  shallow  penetration  theory’’""  is  avoided  in  this  consideration. 
Instead,  it  is  assumed  on  balance  that  the  onset  of  resistance  to  pene¬ 
tration  in  snow  takes  place  at  an  indentation  depth  corresponding  to  the 
centroid  of  the  nose  hemisphere.  That  is,  when  q  =  5D/16.  Under  these 
conditions,  Eq.  Ill -36  becomes 


In 


PS  *  ?  W 

ps  +  f  V«/o 


B  D 


2  *S 


l, 


3  M  D 

T  +  B  r  0g 


(q  -  5D, 16) 


(III -37) 


l  3 


The  first  phase  of  penetration  is  completed  when  the  projectile 
encounters  the  snow-ice  interface.  Here, 


q  *  q°s  ■  n  ~ 

Substitution  of  q  =  R  into  Eq.  III-37  determines  the  projectile  velocity 

q„  ,  which  characterizes  initiation  of  the  second  penetration  phase. 
us 
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In  the  second  phase  of  penetration,  the  projectile  encounters  an 
ice  cover  which  is  assumed  in  the  present  problem  to  behave  as  an  elastic-* 
plastic  solid  in  shear,  and  incompressible  under  hydrostatic  stress. 
Attention  is  directed  to  Figs.  5  and  6.  Using  expressions  which  have 
been  developed  in  Ref.  3,  the  radial  normal  stress  at  the  snow-ice 
interface  is  given  by* 


*  y«p(bS  +  24,sXr  -  0 
- j (r - 0  (III-38) 


Again,  from  previous  work 


=  Of 


r=R 


b2i> 

R2 


fj(t) 


where 


fj(t) 


GT  b2b 


(III-39) 

(III-40) 


*  The 
and 


subscript,  I, 
E  for  reasons 


has  been  omitted  from  the  quantities  Y,  p  ,  p# 

o  * 

of  mathematical  conv<  lienee. 
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FIGURE  5  OYNAMIC  CAVITY  EXPANSION  PROBLEM  FOR  AN 
ELASTIC -PLASTIC,  INCOMPRESSIBLE  MATERIAL 


X  DENOTES  CONDITION  OF  MATERIAL  IN  "EG'0N  'l 
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FIGURE  6  IDEALIZED  STRESS  •  STRAIN  CUHVrS  FOR 

AN  ELASTIC  PLASTIC,  INCOMPRESSIBLE  MATERIAL 
WITH  LINEAR  STRAIN  HARDENING  AFTER  YIELD 
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In  the  snow  region,  the  general  solution  to  the  problem  can  be  written 
as 


(in-41) 


v 


s 


f  (t) 
s 


(III-42) 


Boundary  conditions  are  given  by: 


r  =  a(t) 


(III-43) 


r  =  R(t) 


Pj»  (R  -  v  ) 
s 


(III-44) 


-  Pi  (Vj  '  R)<Vj  - 


v  ) 
s 


(III-45) 


The  quantivies  arj  and  v  at  r  =  R  are  given  by  Eqs.  III-38,  III-39. 

Us^  of  Eq.  Ill— 13  enables  the  quantity,  g  (t),  in  Eq.  III-41  to  be 

s 

expressed  in  terms  of  p(t) . 
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-y*s6jasv,i 


Then,  the  new  relationship  for  a  is  given 


by 


%  *  "  P(t) 


/r 

-L  E  41ri£l  2  _  \ 

r,  E,s  elsj 


dr' 


+  V.  ;  °t/;  (*  -  £) 


(III-46) 


From  the  boundary  condition  of  Eq.  Ill -44, 


/  !fp\.  p< 

r=R  \  PJR+  " 


P*s  1 


(III-47) 


Then, 


f  (t)  = 

S 


f1  -  r* V«  +  f  (t> 

\  ‘V  pis  1 


(III-48) 


Substitution  into  the  boundary  condition  given  by  Eq.  in-45  yields 


o  -a 
rs  rj 


M. =  ■*»(* '  “p  $)  k  7  -  v 


R  -  (l-Q',  )a 


[s  P 


b2b\ 

R  V 


(111-49) 


or 


c„  -  a 
r  s  rj 


r=R 


!  %('  ai-sR2  +  Rafter, 


b'bR 


«  k2 


P  \  *8 


s  P  RC  R4 


(III -50) 


.1 

1 
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Substitution  of  Eqs.  Ill -46,  III -38  in  Eq.  Ill -50  gives 


p(t)  +  2Y  In  -  +  2 
s  a 


T  J ^ 

+  Pig  aAg(2RR®  +  I^S)  +  7“^  Of  (2bb8  +  bs 


SI 

+  c  (2bb8  +  baB)  -  A 
Pjj  P  J\a  R/ 


-\H  f(i  .  +  X  b*bf  (i  -  4} 

2  *l\  pi,/  pl,  "  J  '• 

+  p  t  o!  p  ?  +  #  p.  (abb3  +  bsfi)  (~  -  £) 

I  p  o  p  *p  \R  bJ 


=  D t  V  aL  ^  *  2« »  «  "  O8®* 

*p  \  *s  *s  p  R2  P 


(III-51) 


where 


I  Y  m  fl  -  e-36  ^  ,  £  «•« 

3  \  Pi  /  27  t 


— 1  ' 


■30  Po 


14  -•*) 


(III-52) 
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Solving  Eq. 
surface  and 


HI-51  for  the  dynamic  pressure  variation  at  the  cavity 
subsequent  rearrangement  yields 


p(t) 


2y 


ln»*2/  7e..(2 -f 


+  +  ^  (f  -  j") 

+  w2b6S  +  b£H)  (i  -  £) 


~  2  °*i 


R26  +  7^  <*  b8b 


pi.  P 


Is 

(*- 


+  “  oPpj  b4b2(-~-  -  -L')  +  p 
2  P  *P  \R4  b4  /  i 


+  Vo''  +  ««.%* 


9n  b8  b  . 

20  Wp  ~W  R 


+  e 


b4b2 


P  VS  *P  tf* 


(III-53) 


“•  ,0U°*1W  *»  besn  „erlved  ln  Rsf.  3  „„„  1Uted 

for  convenience. 


2 


27 


^Et  -  -  Et  7] 


(III-54) 


2Y  In  - 
s  a 


2 

3 


Y 

s 


ln 


(III  — 55 ) 
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where 
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Using  Eqs.  III-62,  III-63 


or 


(III-64) 


b 


rS  >ij 


(III-65) 


Now  p a  >  p  ,  and  for  deep  penetration  depths  of  the  order  of  several 
p  o 

projectile  diameters  in  ice,  R3  »  R3 .  Thus,  Eq.  Ill -65  becomes 


^  _  -1/3 

b  =  Ra 

p 


(III-66) 


where 


a  =  (l  -  e 


-38  P  o 


(III-67) 


It  follows  readily  that 


b  =  a  1  -  e 


r1/3  /  .3S  »„ 


o  V  '  / 

-  i- 

PV  \ 


.  -1/3 


(III-68) 


or 


-  -1/3  _  -1/3 
b  =  aa  a 
s  p 


(III-69) 


where 


.  L 


-36  Do.\ 


s  s 

H 


; 


(III-70) 
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Introduction  of  these  relationships  in  Eq.  III-53  yields 


2..  /.  '3PsP°s\ 


p(t)  =  -  —  Y_  In  ^1  -  e  s  ——  )+  ~  r^E 


3  s 


Os  \  2_ 

nj +  27 "  "ts 


2  % 


2  • 
a  a 


a  2  •• 
(a  a 

+  2aa2) 

(l- 

i) 

s 

\a 

R  / 

(a  )-1 

(a2  a  + 

2aa2) 

(i- 

-A 

s 

\a 

b/ 

-  Q 

f  (5  >** 

(a  )-1 

2  * 
a  a 

2 

(4 

p*s 

p  p 

s 

\a4 

*) 


+  P  +  a  P  (5  )"2/3(c i  )-2/3a2  +  P£  <*4  W  )  ^  a2 
I  p  o  p  s  *p  *s  s 


-  2 Pi  or*  a?  (5  ) 

*p  p  p 


(S  ) 
s 


+ 


H/fi 


/-  x-3/3  >2 
(ff  )  (a  )  a 
s  p  s 


(III-71) 


Now  let 


P 

I 


+ 


(III-72) 
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Then 


33 


It  has  been  shown  in  Ref.  3  that  an  expression  for  dynamic 
pressure  at  the  spherical  cavity  surface,  such  as  the  one  given  by 
Eq.  III-75,  can  be  employed  to  derive  an  equation  of  motion  for  pro¬ 
jectile  penetration.  In  the  resent  case,  the  following  relationship 
is  obtained 


( III -78) 


31 


The  usual  process  of  Integration  yields 


In 


+  2/3  C  q2 
+  2/3  C  v2 

2  °S 


”0,) 

3  M/A  +  C  D/3 


Solving  for  the  terminal  penetration  depth 


(III-79) 


=  <lr 


3 

+  - 


M/A  + 


C  D/3 

l 


In 


II 


+  2/3 


C  vs 

S  °S 


(III-80) 


II 


Here,  q_  and  v„  refer  to  the  penetration  depth  and  exit  velocitv, 

°s  °s 

respectively,  characterizing  projectile  transit  through  the  overlying 
snow  cover.  These  quantities  are  obtained  from  solution  of  the  first 
phase  penetration  problem. 
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Appendix  A 

LISTING  FOR  LARGE  DEFORMATION  THEORY  COMPUTER  CODE 
THAT  INCORPORATES  STRESS  REDISTRIBUTIONS  DUE  TO  FRACTURE 


Preceding  page  blank 
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JOB  485 


//LIFTER  JUH  '  A5 19*** ,602  » 1*9*  * • SIOHU* 

//•SERVICE  tIST 
//•POINT  COPIES*? 

//$rpPl  EXEC  FORTHCIG 
//FU«T.SYSIN  CC  * 

C  »**t ••«•**•*****«********««**»****#«**********»*•««*«•*•«**«**** 

c 

C  VFRSICN  ESPI....TWO  OfMENSIONAt  AXI SYMMFTTRIC  ‘CANDIA*  COD’:, 
f. 

C  THIS  VFRSICN  COMPAYIRLF  UN  ANY  IBM  160  MACHINE 
C 

C  LAnG'»AGF....FCRTRAN  (IVI  h 

c 


c 

C  *.  «>  MM«»MM******»»*«**n****M**t********* 
IMPLICIT  REAL*8  <A-H,0-Z,»),  INTFGFR  (t-N) 


COMMON 

/CCMl/ 

AR  ( ?4  ,46 ) * 

ARH 

(24,46), 

ARDMI 24,46 ) , 

1 

A/  (24,46), 

AZH 

(24,46), 

AKIHI  24,46), 

2 

A A  ( 24  ,46 )  , 

AV 

(24,46)  , 

AP  (24,46), 

i 

ATOP  ( 2  «,4o) , 

ATZK  24,46)  , 

ATRZ ( 24,46) , 

4 

ATTT ( ?4  ,46)  , 

S 

A<RP(24,46)  , 

AS? Z (24,46)  , 

ASRZ ( ?4,46 ) , 

ft 

A5TT( 24,46) 

ciimon 

/CCP?/ 

HR  (46  ,  6)  , 

BRH 

(46 ,  6 )  , 

BP OH (46,  6), 

1 

R /  (4o,  6), 

SZH 

(  46 ,  6)  , 

RZDM(46,  6), 

A 

t 

HA  (46,  61, 

HV 

(46,  6), 

RP ( 46 ,  6), 

) 

BTRPUft,  6)  , 

RTZZ(46,  6), 

BTRZ (46 ,  6), 

4 

BTTT ( 46 ,  6), 

6 

HSRP 1 46  ,  6)  , 

BSZZ<46»  e.)  , 

ftSRZ  (46,  6), 

r 

ft STT (  4ft ,  ft) 

-  -  . 

COMMON 

/CC*T/ 

TR  C’4,46), 

TRM 

( ?4 ,46  )  , 

TCOH( ?*,«,4) , 

l 

T/  (?4,46), 

TZm 

(24,46)  , 

T7ni((24,46), 

1  S«(4ol,  SZI46),  TSR  (  46 1  »  TSZI46),  ZM(2,46),  SPTM(*5), 

4  All?),  A2I?),  A3  12  )  *  AAC»),  At  EN(  ?  I «  AMAT 12 ) »  A  10(2),  4NOI2), 

•>  RK( 2 )  *  CAPF(?»,  ff>UT(?),  OR  I  2 )  *  D? ( 2 )  »  RAD(2),  RHI'i  (21*  *»D0T(2I, 
(■  RHfV.(2),  PMC  T  (  ?  ) »  Thpu(2>»  VT  C  2  > ,  ZDi)H2>,  ZV  f  2)  »  VCJNt  ?  ) » 

7  uQ»DTN,C  lP.f>TH»s,!HMIN,ZnlVTf  TVMF  »ChFK0»F4CTP 
C  JM  MON  /ICC*!/  I  PT  A  <  24*461 

ClmmoN  /ICCM2/  l  PTf)  (46 »  61 

COMMON  /ICCM3/  I SHP ( ?  I »  JM(?)t  KM(2>,  KM|P, 

I  I P'JNt I PTI » ICT>,  ICFNT,  !MAX,IPRNT,ICYCL»  IMAP, NTT, K INT ,KMN, KMX 

r  ***** 

r 

r  SFT  ARTIflClAl  VISCOSITY  COEFFI C  IF NT  VAUJF 

c 

C  ft*  HA* 

C0=C.1 

HEAD  (*,!(•')  10FN  T 

loo  r via ma r  (  no) 

c  >»*♦* 

c 

C  10 ENT  *  l  NFR  RIJN  Ml  SAVF 

C  *  2  Nf K  RUN  SA VC 

f  -  ’  PfSUMf  NO  SAV* 

r  s  4  cf^upf  save 

r 

C  *«-*»* 

!c  (  IOFM.Cc.T!  CM  L  »SM? 

IF  ( I0ENT.CF.1 I  oO  TO  -P? 
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c  ***** 

c 

C  READ  IN  CATA  FRCM  CAROS 
C 

C  ***** 

CAU  INPTT 
C  ***** 

C 

C  DU  ALL  INITIAL  MODIFICATION  OF  DATA.. .UNITS  FTC 
C 

C  **  **« 

CALL  INITT 
f  ***** 

C 

C  CREATE  TFE  FINITE  DIFFERENCE  GRID  WORK  AND  ALL  DISC  FILFS 

C 

C  ***** 

call  gricc 

c  ***** 

c 

C  WRITE  OUT  THE  SALIFNT  ASPECTS  OF  THE  PROBLEM 
C 

C  *»  *** 

3"?  CALL  INWRR 
C  ***** 

c 

C  CYCLE  THROUGH  THE  RUDY  IS  INITIATED 
C 

C  ***** 

CALL  LETGC 
C  ***** 

c 

C  DO  ALL  NECESSARY  COMPUTATION  FOR  THf  TARGET 
r 

C  ***** 

CALL  viUTT  (  l,IPTAf  AP , AZ, ARH, AZH , ARDH , AZPH,  ATRR ,  AT?  Z ,  ATR  7  , 

1  ATTT,ASRR.ASZZ.ASRZ»ASTT,AA,  AV.  AP,KMm,KMN»K INT, 

2  JMX) 

C  ***** 

c 

C  BEFORE  GCING  TO  THF  PROJCCTUF  00  All  PREPARATIONS  NECESSARY 

C 

C  ***** 

JJJ*J  M ! 2 ) 

00  304 

BRHll  , J)  «  ARH1KINT, J! 

BZ  HI  1  *  J  I  «  AZH 1 K INT  *  J  ) 

BROHI 1  *  Jl  *  APOMK  INT  *  J  ) 

3r4  BZOHl  1  tJ)  *  A  ZDH  ( •<  I  NT  *  J) 

C  ***** 

C 

C  NOW  00  AIL  THE  COMPUTATIONS  FOR  THE  PROJECT  ILF 

r. 

c  ***** 

K2  *  K«(2) 

CAU  GUTT  (  2,  IPTB.BP  ,  BZ  ,  8RH,  BZH,  BR  DH , «  ZOH,  BTF.R  ,  T  ?  2  ,RT<>7  , 

1  HTTT,:>SRR,BSZZ.BSRZ.BSTT,3A,B'',BP,*2,  2,*M*t 

2  JM ( ?  )  ) 

C  ***** 

C 
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r  OtfC<  If  IT  IS  TIMf  TO  v»-ilT 
r 

i  »»•*** 


If  I  ICvCl-lMAX)  r",K',T"! 

*>!  If  MI*'fM.fC.?».fl».l  !C»;f'T.=0.4n  CAl  L  SSAV 

rft'ips 

S'P.M)iiT|K*  PSCf 

i^nicn  r€al*4  < A-mn-; . *) ,  inteofr  u-ni 


Ctiv»»s;r.  /rc»l/  A«  (24,46), 

l  A/  (  24  ,-*•> )  « 

>  4A  I  P4 ,46 I < 

*  ATfO  <  >4  ,4<>)  , 

4  AT  T  T  (  N*<o|  , 

6  ASP 2(24, 46), 

*>  AST  T(  ?+  «4f> ) 


APH  12 4,46), 
4/H  ( ?4  »4:i  I  « 
AW  12 4,441, 
AT?Z<24,46) , 


AR4M<  2  .,46) , 
A20‘M^4,4M, 
AP  IM.4M, 
Wall  >4,44)  , 


AS/ / ( 24*  46 ) ,  AjPM24,44), 


CO  A  A.j‘4 

/CCM7/  fl»  <46, 

ni  ♦ 

RRrt  <46, 

6)  , 

HP OH <46, 

6) 

1 

R7  (4*1, 

)  * 

R7M  <46, 

6), 

H70MI 46. 

4) 

1 

t 

HA  <  66 , 

ft)  ♦ 

RV  <46. 

6 )  * 

P  °  <  4  n ,  6) 

, 

RTRR (44 , 

M  t 

«T»7<46, 

6), 

R TP ? <46, 

6) 

** 

RTTT  <  46  , 

M  t 

. 

i 

BSRO  <  44 , 

i>l  t 

RS7ZI46, 

6)  , 

RSR  7 ( 46 , 

6) 

UlMRON  /ffMV  TR 


BSTT  ( 4t>  *  6) 


T/ 


< ’4 , 4o ) , 
< 24,46), 


TPM  < ->4,461, 
f  ZH  I  74*46 ) « 


TP  'Vtl  24, 4o  >  * 

T/nm?4,46)  , 

Sw<46>.  S7<46),  TSPI46),  TS7<46),  ZM<?,46),  SPTMIS), 

Al ( 2  I »  A2I2J,  4J(2),  A4<2),  AIFM/J,  AMATI2),  AH'U<>),  A‘itj<2), 

■  R<<?1,  C  APP 121,  ( f'cT(2)  *  n  K  ( .? ) ,  ,)Z<2),  pad!?),  phu<2).  ROOM?), 

i  RMllu(?),  RNC!T<?)t  TWM'jMI,  VT<?J,  ZOOTI?),  ZVC>>,  VCiJM’), 

CO,MTN,rTh,nTHW,rTM|N, hunt, tvmf.chfkd.factp 
C.JMON  nc.c»l/  |PTA(  24,46) 

COMiUM  /ICCM2/  IPTHI46,  6) 

COHJN  /ICCKR/  I  SH°  ( 2  )  *  JMI?),  KM{2),  <M)0, 

IR'JM*  ITT  If  IDTPtlOFNTtlMAW* lPKNT.tr  vet , IMAP, NT T , KINT ,KMN, KPT 
RF  T'J«\ 
f  NO 

Sopp.KIT  INF  SSAV 

implicit  pfai*h  < a-h,oz, n ,  jntofp  c i-m » 


r  /r cm i /  ar  <24 

1  A/  <24 

2  AA  <’4 

1  A  f  S  K  <  9  4 

4  ATT! <24 

4  ASF  PI  24 

'■  ASTT  (  24 

CO i*'  AUN  /CCP?/  BO  <46 


46  )  , 
44)  , 
46)  , 
46)  * 
46)  , 
46 )  * 
46) 
61  . 
6 )  , 
6  )  * 
4)  * 
4)  , 
6)  , 
6) 

46)  , 
46)  , 


4R»t  124,46), 
AZH  !2 4,46), 
AV  <24,46), 
AT/Zt?4,46) , 


AP»M(?4,4o), 
A/DHI 24, 46 ) , 
AP  <24,46), 
ATR  7  ( 24,46 ) , 


A  SZ7 ( 24  ,46 )  ,  ASP  I ( 24 ,46 ) , 


HPH  146,  6), 
RZH  (4e*  6) , 
11V  (46,  6), 

BTZ7 146,  6) , 


RRDH<46,  6), 
RZr)H(46,  6), 
RPI46,  6), 
RTRZ ( 46  *  6), 


R?  <46, 

■>  RA  (44 

>  i)TFR»46, 

4  RTTTI 46 i 

4  tt  SKIM  44 

6  P  ST  T  (  4fc  i 

CUMMCN  /CC*3/  TO  <24,461,  TOH  <24,46),  TR0H<24,46), 

1  T  7  (24  ,46),  TZH  <24,46),  TZ0H<?4,44), 

2  SM46),  S  2  ( 46 ) ,  TSRI46),  TS/146),  ZM(?,46),  SPTRtR), 

4  A l ( 2  ) »  A2( 2 ) ,  A  3(2),  A4 ( 2 ) ,  A|fN<7).  A*  AT ( 2 ) »  AMOI 2  I ,  ANU<2), 

4  R*<2),  CAPPI2),  OOTI’I,  HR  <  2) ,  JZ ( 2)  ,  R  AD  t  2  » ,  RM0I7),  RDOTIZ) 
4  MUG  12),  RHCT  (  2 )  *  TWP'H.?),  VT  <  2  ) ,  Zr>t)T(2),  ZV<2),  VC0NI2), 

!  r  J,OTN,rTH,DTHV,,PTMIN,Z')INT,TVHF,r.HFKn,FACTR 
C  lMMCN  /SfCPl/  I PT  A  I ?4 ,46 ) 


RS  Z  7 ( 44  »  6),  8SR7  <46,  6), 


NOT  REPRODUCIBLE 
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,~ir~>oor)nnr->Of_>o  o  r> 


COMMON  /(COM?/  IPTB(46,  61 

CCNNON  /ICCM3/  ISHP(2),  JN(2),  KM(2),  KM(Ot 
l  I  RUN* ICT l» IDT2  tIOENT *INAX»  J  PRNT, ICYCL » IHAP*NTT  ,KINT ,KMN, KMX 
RETURN 
END 

SUBROUTINE  INPTT 

IMPLICIT  REAL *6  (A-H,0-Z,S),  INTEGER  (l-NI 


CUMMUN  /CCM 1/ 


1 

2 

3 

4 

5 
4 


AR  (24,46), 
AZ  ( 24,461 , 
AA  ( 24*46 I , 
A TRR( 24  *46 )  * 
ATTT(24,46), 
ASPP( 24  ,46)  , 
ASTT ( 24,46) 


ARH  (24,46), 
A2H  (24,46), 
AV  (24,46), 
ATZ  Z  ( 24,46) , 


AROM( 24,46) 
AZDM( 24,46) 
AP  (24,46) 
ATRZ( 24,46 ) 


ASZZ(24,46>,  ASR7 ( 24,46) 


CCMMON 

/CCM2/ 

AR  146,  6), 

BRH  (46, 

6), 

BP DH(  46,  6) 

1 

RZ  (46,  6) , 

BZH  (46, 

6)  , 

BZDH(46,  6) 

•> 

BA  ( 46  ,  6 ) , 

8V  (46, 

6), 

BP { 46,  4), 

* 

RTRRI46,  6) , 

BTZZ(46, 

6) , 

BTRZ ( 46,  6) 

4 

BTTT ( 46  ,  6) , 

5 

BSPR(46,  6), 

BSZZ(46, 

6 ) , 

BSRZ ( 46,  6) 

6 

ASTT (46 ,  6) 

CUM  MUM 

/CCM3/ 

TR  (24, 4o), 

TRH  (24, 

46)  , 

TROH( 24,46) 

1  TZ  (24,46),  TZH  (24,44), 

2  SR ( 46 )  •  SZ ( 46 ) »  TSR ( 46 ) ,  TSZ (46) , 

4  At  (2 ) ,  A2(2t,  AJI2),  44(2),  Al.EN(2>,  AMATI2), 
r»  BK(  2)  ,  CApF  (  2 )  ,  Cr»nT(2),  DR(2),  DZ(2),  PAD(2), 


TZDH(24,46), 

ZM(? ,46) ,  SPTM(5>, 

AMD(2),  ANIH2), 

R  HD  ( 2 )  ,  ROOT  ( 2  ) , 


4  RHOG(2>,  RHOT  (  2  )  ,  TWPUI  2  ) ,  VT(2),  Zr*OT(2),  ZV(2),  VCON(2>, 

7  C0,DTN,CT6,DTM4,0TMIN,Z0IiMT,TYMF  ,CHEKD»F  ACTR 
COMMON  /  ICCM ( /  1PTA (  24, 46) 

CCMMCN  / 1 CCM2/  !PTB(46,  6) 

COMMON  /  ICCM3/  ISHP(2),  JM(2),  KM{?),  KMJO, 

1  IR'JN,  ICTl,  (DT2,  l TENT  ,  IMAX,  I PRNT ,  If.  YCL  ,  IMAP ,N  TT,K  IN T  ,XMN ,  KMX 
"tCR*5 
C  */ **-«* 

C 

r  P  F  AO  IDENTIFIERS: 

c 

c 


IRON  :  RUN  I.O.  NO 

IOT1  :  MCNTH-OAY  IFG,  1103  -  3RD  oOV ) 
I0T2  :  YEAR  (EG.  60  -  PQR  I960) 

FAC TR  :  FRACTION  CF  DIAGONAL  USEf  ru  FIND 


THE  TIMr  STEP 


**  **** 
READ 

****** 


(NCR, ICO)  IRt»N,!DTl,I0T2,FACTR 


RFAO  CuNTPLl  VARI ARIFS 


IMAX  -  MAXIMUM  NO,  OF 
I PRNT-  MAXIMUM  NO.  uF 
IMAP  -  MAXIMUM  NO.  OF 
NTT  -  NT.  CF  SPECIAL 
SPTMINTT)  SPECIAL  PRINT  TIMES  IN 


CYCI.FS  BEFORE  PROBLEM  TERMINATION 
CYCLES  BETWEEN  PRINTOUTS 
CYCLES  BETWEEN  GRAPHICAL  OUTPUTS 
PRINT  TIMES 

MICRO-SECOND* 


****** 

READ  (NCR,  IT’D  IMAX, IPRNT, IMAP, NTT 
RE  AD  (NC», 102)  (SPTM(K),  K« l, NT  T  ) 

♦«  **** 

MAO  MATERIAL  PROPERTIES 
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r 

c. 

r 

C 

c. 

c 

c 

r. 

C 


r 

f 

r. 

r 

c. 

c 

r 

c 

r 

r 

r 

c 

c 

r 

r 

C 


r 

r 

r 

r 

r 


AMAT(?)  -  MATER  I  At  NAM*  (4  CHARACTERS  MAXIMUM) 

AN'J(2)  -  PCISSCN'S  PATIO 

CAP£(2)  -  YOUNG'S  XOOtILUS  -  PSI 

RHO(2>  -  CENSITY  -  GM/CC 

VT(7)  -  TENSILE  Yfnn  STRENGTH  -  »SI 

Al-A4<2>-  SPARES  ■  f 

******  * 

on  *oo  j»i»2  '  , 

READ  (KCRtlPS)  AMAT(J),ANU( JI»CAPE(J».RH(HJJ.YTIJ) 
RE  AO  (NCR, 1*4)  A1  fj  >  ,  A?U  J ,  AT  ( J  )  ,A4(  J) 

****** 

E EAO  GEOMETRIC  PPOPERTfES 
t  SHP( 2  )  -  SHAPE  FACTOR 

BODY  l  (TARGET)  >  V  ,ROO 

2  .PLATE 


RuOV  2  (PRLJFCTUE)  *L  ,PCO 

7  .PLATE 


ROOT (2) 

-  INITIAL  VFLCCITY  IN  R  DIRECTION  -  FT/SEC 

ZOOT ( 2 ) 

-  INITIAL  VELOCITY  IN  7  DIRECTION  -  FVSEC 

-  no*  oe  cells  Yin  r  direction 

JM(2) 

KM  (?) 

-  ISL.  OF  CELLS  ?H  z  OIREfTIUN 

OR  (  ? ) 

-  SIZr  (if  CcL LS  IN  R  OIPFCTI0N  -  INCH 

ozm 

-  SIZE  r E  CELLS  IN  Z.  DIRECTION  -  INCH 

*  * 

*«*¥■ 

HOO 

RE  AO  (NCP.10S)  I  SUP  (  J  )  ,  ROOT  (J),ZOflT(J),JM(J),KM(J),DR(J),nZ(J) 

1 

FijRMAT 

(31 10, FI". S) 

!■". 

FJRMAT 

(AH'*) 

1*7 

format 

(EE1C.0) 

n 

FORMAT 

(6 X,A4,4E10.F) 

104 

FORMAT 

(4E10.0) 

r  s 

r. JR  MAT 

RET  IRN 

HNO 

( lie, 2EK. -7,21  n,2Fl'’.') 

SOMKJtlTIAE  IMTT 

*  <X  *  t  * 


THS  SOMRctmE  PREPARES  THE  (mPUT  DATA  FuR  PRINTING  ANO  SUBSEQUENT 

iKp 


*<  *  *  « 

IMPLICIT  RR4l*P  (A-H.O-Z.t) 
CLM.MUN  /CLP 1/  AM  (24,46), 

1  4/  ( ?6  ,44 (  , 

7  A  A  (74,46), 

)  4TPM(?4,46), 

.  A  T  T  T ( 74  ,46 )  , 

6  ASF  R  (  24 ,46  )  . 

<•  4  ST  T  (  24  ,46  ) 


,  INTEGER  (l-N) 
ARH  (24,44), 
AZH  (24,46), 

AV  (74,46), 
ATZZ(?4,46), 


AROH( 24,46) , 
4Z OH (?4, 46), 
AP  (74,46), 
ATP  ?  (  ?4, 44) , 


A  SZ  Z ( 24 ,46 ) ,  AS0,( 24,44)  , 


COMMON 

/CCM2/ 

l?R  (  46  ,  4  )  , 

oph  (46, 

6)  , 

ROOH( 46,  6), 

1 

RZ  (46  ,  A)  , 

HZH  (46, 

6 )  , 

BZ  OH ( 46 ,  4), 

7  A  '4a,  6)  , 

nv  (46, 

6)  , 

RP(46,  6) , 

■X 

MTRR (44  ,  6 )  , 

IT  Z  Z ( 4  3 , 

41. 

RTRZ(46,  6), 

4 

0  T  T  T ( 4  4  ,  4)  , 

4 

RSRR(46,  6)  , 

BSZZI46, 

6)  , 

P SR Z ( 46 *  6), 

t 

«STT(46  ,  6) 

CUN  MON 

/CC“V 

TP  (74,46), 

TOO  (?4, 

46)  , 

TR  OH ( 24 , 46  )  , 

l 

T  7  (’4,46 ), 

TZH  (74, 

46)  , 

TZOH ( 74,46) , 
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2  SR146),  S2  ( ^6>  TSPC40),  TSZ(46)/,  ZM(?,46>,.  SPTJH5), 

4  All'll  A2(2J,  AT(2),  AM  21,  U  ENI  2  > ,.  A«A  T(  21,  AMU<  2.1-,,  ANU(  ?1, 

5  8K(  2 )  *  CAPE( 2  I  »  0)01(2).  OR ( 2) ,  OZ (2 1,  RA0C21,  RHO  { 2 )  ,  ROOT  I  21. 

6  RHOG(?),  RNOT<?);  TWMU<21,  VT(2j,  ZOOT  (2 » .  ZM2(,  VCON{ij, 

7  CQ.OTN  ».CTE»  CTHh,f,TMI  N.ZO.INT  ,  TYME  »CHFKO,F  AGTR- 
COMMON  /1C0M17  iPTAI 

CCMMON  Al COM2 /  IPT8(46,  6) 

C-CMMON  /IC.CM3/  I  SH°  (2  ) .  JM(2!  .  KM  (21,  KMIU, 

1  IRON,  I CT  l,  IDT 2 »  IDFNT  ,  IMA  X,  I  PRNT ,  IC  *CL  »  IMAP  iNT-T.,  K INT-, KMN,  KMX 
no  301  J  *1.2 
RHUuU)-  RFC! f  J  ) 
p.H(HJ)  =  9.354320-05  *  RHO(.J) 

C  *.'*>(* 

C 

r:\in,  rhoo  is  in  (c.m/cci 

r  PHu  .IS  IN  1LB-SFC2/IN41 

r 

fi,  **  *«■'* 

AMU ( J )  =  (CAPFIJ  )*c,5)  /( 1.0+ AMJUn 
BM'JI  =  CAPE(-J1/H.0*(  l. 0-2.0* AMD (J  )  )  ) 

COOT!  J)=nSOP.T(  (P«(  J),+  4.0*AMii(  JJ/i.^J/RHOU)  1*1. 00-26 
VCON (  J  )  =  COOTU  1/1 .20-05 
C  *rr** 

c 

G  UMTS:  AMO.fiK  ARE  IN  (PS  I) 

C  CDCT  15  IN'  (INCHES  pfp  MICRDSCn 

G  VCCN  IS  IN  (FEET  PER  SEC) 

r 

C  *t.  tt  n. 

AlENIJl  =  CZ(J>  *  KM ( J ) 

R.W(J)  =  CPU)  *  JM(.J) 

■f,  *Rr.V«- 

C  NCvf  CHAWuE  JM  AND  KH  TO  STANO  FOR  NO  OF  t.PIO  RONS  A  NO  COLIMNS 
C  INSTEAO  CF  NU  OF  -i<MO  CPUS 

C  ♦*■»<•* 

JM  J )  =  JM(j):  +  l 
Kit  J)  =  KHI  J)  f  1 

/.v(j)  =  zncrrji 

zivii(j)  =  1 2.00 -06  *  ccorui: 

c  ***** 

C  UNITS :  ZV  IS  IN  (FT/srci 
C  ZCT  T  IS  IJ  '(luC'l  PFR  MICRCSFC) 

c  *<■*** 
c 
c 

C  NFXT  COMPUTE  VAI UES  FOR  TIMF  ST^P  IN  MICROSEC 
r 

OTHvJ  =FACTR*OSORT(i)R(  J,)**2  +  r7(J)**2.L/  COUTH) 

C  *A*R» 

IF  (J-l)  3C0, 300,201 
Vf"  OIH=OTHVi 
CONTINUE 

r  *  r  .*  Hi « 

C  <jC*  SELECT  TEF  SMAtlFP  DF  OTM  f.  f)THV»  TO  BP  flip  1|“F  ST'l-o 
C  ***!,.« 

IF  <  OTI — CTHw )  >''3,303,302 

?r  2  or h=9Thvs 

303  DT MIN  = 'CTH  *  0.1 
KI  ivT  =  K'N  il) 

KM  1 0  =  K!i\T  +  J 
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.'V  o  o 


C—  COMPOTE  TK  INTERFACE  VELOCITY 

f  *«9Ci; 


dint*  izooti  u*iiHi'iii*r.wmi  ♦  Z(wi2I*pmoi2i*c  n'Tmi/iwHOM  i 

I  CCCT  (  l  I  «•  RMU(7)*CO(!T|2M 

P-  TUP  N 
FNO 

•;oa»j'jT  iae  inwrr 

IMPLICIT  RFAl«R  U-O, *,-?.♦>.  IMFGr.fi  { I -M 
COMMON  /CP**!/  AR  I  ->4,46),  4PM  (?»,4o) , 

AZH  <2 *,44), 

AV  (24,46), 

AT77(?4,46), 


1  47  ( ?4,46 t. 

»  M  I74|44li 

■S  ATF»(  ’4,46), 

4  AtTT  ( ?4  ,44  I  , 

1  A<RR(24,46), 

6  AST  T(  74, 46) 

COP  MON  /(T*?/  MR  (44,  Dll 
l  «7  (46,  6), 

/  4  4  (46  ,  6)  , 

•  HTUH4(,,  o>, 

4  PTT  T  ( 4»»  ,  6)  , 

w  Q  SPP.  ( 4o  ,  41  , 

4  RSTT ( 46 ,  6) 

r.JlirtUN  /CLF3/  TP  (24,44), 
T  7  (74,4b), 


AROrll  ^4 , 4<i  )  f 
*710(74,46), 
AP  C»4,44), 
A  TP  * ( 74,44 )  , 


ASZ  7 1 ?4  ,46)  ,  ASR7J >*,4o) , 


BRH  (46,  b), 
AZH  (46,  6), 
(V  (46,  6), 
A T7 7. ! 44 ,  6), 


BP  Did  •»{,,  6) , 
P.7  OH  (  44,  6), 
PP(46,  (), 
RT»’<44,  4), 


)S7M4o,  6),  HSR7(',6,  4), 


TRM  (24,46), 
T»(1  (24,46), 


TP 1H( 74,46), 
T70M(.»4,46)  , 

2  SR  ( 46  )  ,  S2  (  46  )  ,  T  SR  <  44  )  ,  TS7  (46)  ,  2M(?,44),  SPT.M(E), 

«  Al  ( 2 ) ,  A2(2),  A3I7I  *  A4( 7) »  ALEM7 ) ,  AMAT(P),  6*1(2),  ANO(’), 

6  DM2),  C  APE  ( 2 ) ,  OCTI2I,  PPI2),  07(2),  «A0(2I,  RHOU),  P')0T(2), 
*.  PHJP,(?),  PHCT  (2  ) ,  TMP’*(  2  ) ,  VT  (  2  1 ,  20t)T(2),  !V(2l,  VCi)N(’), 

’  Ct3,OTN,r.TP,OTHW,rTM|N,Zl)!NT,TVMF,CMFKO,FACTft 
CCMMCN  /ICrMl/  l 6T  A ( 24,46 ) 

C CM MON  7ICCM2/  (PTH(46,  6) 

COMMON  /ICC.M37  ISt-°<2),  JM(2),  KM(2),  KM(0, 

1  («  J.M,Trm,(OT2,|0FNT  ,IMAX,!PRNT,ICYCL,  IMAP ,  NT  T  .HINT  ,MN,  KMX 
NpR*6 
r  «,«<*« 

f 

C  WRITr  PROBLEM  IDENTIFICATION 

r 

c  ***..* , 

W°  ITfc  ( APR, 7f  r )  IRON,  1 0  T  1  ,  lOT?ji  I  Of  NT 

r  *«*’»** 


UTF  MATERIAL  PRlIPFPTlES 


r  <•■*«** 

WRITE  (NPR, ’PI  )  AM  AT ( 1  ), AMAT( ? ) , RHOGt I » ,PHUG(  »  I  ,RM(,(D  ,  RHId  2) 
WRITE  (APP,?P2)  AN.Jd  )» AN1  d?  ) , CAPE  (l ),  CAPE!  21  »AHi)(  1 1  ,AMU(  2) 
WRITE  (f  PR,  2C  3)  BK(  1  )  ,f)K(2)  ,C00T(1)  ,C nnT ( 2 ) , VCPN !1  )  ,VfON(?) 
wKITL  (NPR, 204)  VT ( l ) , VT ( 2) , M (  1 1 , A l( 2 1 

C  ■nlttl 

f 

r  i.rf  Mr  TRI  C  PROPERTIES 


r  t « ***<- 

..RITE  (NPR, 204)  ALT  M ' ) »  ALCW( 2 ) »  R  AO  1 1 ) »  RAO  ( .? ) ,  OR  ( l )  ,DR  (  2 ) 

WRITE  (  N'PR  *2Cfc )  n?(!).0,.»7),KM(l)fKM(?),JM(l),JM(?) 

wr  1  T  F  (NFR,20?)  ROOT  (  1  ) ,  ROOT  ( 7  ) ,  ZP;)T(  !  )  ,  ZOOTI  ? )  ,7.  VI  1 )  ,  Z  V(  2 ) 


,c; 

r; 


CONTROL  047* 


'***♦%* 
WRITE 
WRITE 
**  *** 


FAC:TR  ,  IMAX,  !HAP,IPRNT»NTT 
c  *****  '  JN??*?09»  f  SPTM { KKT ) »KKT= I ,NTT) 

C  ,  ‘ 

C  FORMATS 

C 

C  ***** 

20*  tF0RMAT  TM*.,70X,  •  TWO-CLAGR*KGI  AFT  COMPUTATl  ON*  ,/, 7 IX,  *  RUN.  NO 


H,<  CATE  ^  SI4.S  19»I2,10X,  •RESUME  CODF  *-»  U 


'BODY  B V* / / ) 


,  *  ■  - - 

//// »40X»* UNITS* *?9X9,BndY  AJ«1*X, 


SOOULUS  , IcwIkISx *5s|?  'l»r?niS^l;,3i?,2"0->’,/l5X'??H™"*“!,S  »■ 
S' ;  13* !“sf  y  ?5;»f  {«' f  Jf '  !C20-<"  «•••»«"''  TV  MODULUS  r.MU,  • , 

^FORMAT  (5X,  'BULK  MODULUS  (SKI  18Xt  'PSI  '  ,  17X.2D’0  6  //  «;y 
S  .SOMNC  SPEEC  (CDOT  )S1TX, '',IN/M!CRO-SEC  'j8xf°t0‘6,//,5X’ 

S  2D2C.6, /  »40X  ,  ‘FT /SFC'1 , 

S  14X,2D20.6* /J 

FORMAT  (5X, 'YIELD  STRENGTH  ( YT )  *  ,16X,  •  PS  17X . 2D70  a  / 
format  .  Iwa(. 

FORMAT  (SX,«0Z*,33X,*lNGH',IfcXf2020.6,//,5X,'GRIO  POINTS  l  (KM) ' 

S  37X, 2120  ,//,5X, 'GRID  POINTS  R  (JMI ♦  17*  4,,V ,  '  11  ’ 

s'  ‘CCMRCL  tATAS/f.’J"JC‘6’/,4<',''fr,SEC'’141t'2'’2'’-4'"/'2*‘, 

3  |"Xu.mSc"x?  1  PRI4TS '  >1  Ax, I10//5XTUSPECI AL  PRINT  TIMFS 
209  FORMAT  (5F15.5) 

RETURN 

END 

SUBROUTINE  GRIHO 
IMPLICIT  REAL*8'  ( A-HiC-Z  , i) 

COMMON  /CCMl/  AR  (24,46), 


20  Z 


203 


294 

205 


206 

2C7 


AZ  (24,46) 

AA  (24,46) 
ATPR  ( 24-, 46  ) 

A  T TT <  24  ,4o ) 
ASRR ( 24 ,46 ) 
A:STT(24  ,46) 
COMMON  ./CGM2/  BR  (46,  6) 

1  B l  (  46  ,  6 ) 

2  BA  (46,  6) 

3  BTRR ( 46  ,  6) 

4  8 TTT ( 46  ,  6) 

^  RSRR ( 46  ,  6) 

6  B ST T (46  ,  6) 

COMMON  /CCM3/  TR  (24,46), 

J  TZ  (24  ,46)-, 

2  SR ( 46) ,  SZ(46),  TSR(46), 


INTEGER  (I-N) 
ARH  (24,46), 
AZH  (24  ,46), 

AV  (.24,46),, 
ATZ  Z  ( 24 ,46  )»• 


'AROH(  24,46)  , 
AZDH( 24,46), 
AP  (24,46), 
ATRZ (24,46), 


»'  ASZ  Z  (  24 ,46  ) ,  A'SRZ  (  24,46)  , 


8RH  (4b, 
BZH  (46, 
B.V  (46, 
BTZZI46, 


6) 

6) 

6) 

6), 


t  B  SZ  Z  ( 46  ,  6) 


TRH  (24,46)  , 
TZH  (24,46), 
TSZ ( 46 )  , 


BRDHI44,  6),. 
BZ0H(46,  6), 

SP  (46,  6T, 

B  TR  7-(  46 ,  6 )  , 

BSRZ(46,  6), 

TROHt  24,46), 
TZDH( 24,46 ) , 

ZM( 2,46) ,  SPTM( 5 ) 
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4  41(28,  42(2),  43(?».  A4(2),  ALFN(2),  AMAT  ( 2)  *  AM'I(28,  ANU<?8, 

5  8M2I,  C  APf  { 2 )  »  cnnr(29,  I«Ri2»,  OZ  » 2  »  V  RAO(? ) •  PHl*(?>,  on(>T  C?  » , 

6  RrtOC<2),  PHOT <  2 ) »  Twmm,  Y T(2),  ZD0TI2),  ZVJ2I,  VC0NC>», 

7  LQ,OTNtCTH,OfHl«tnTH!N,ZOTNT,TYMF,CHPKf>»FACTR 
COMMON  /  ICCM1 /  IPTAI 34,46) 

COMMON  /ICCM2/  tPT(5<4A,  6) 

COMMON  / IC0M3/  ISHP  I  ? ) »  JM(2) ,  KM(2),  KMIO, 
l  IKON,  ICI  It  IDT2,  IDFN1  ,I‘1AX,  1PRNT,  ICYCl ,  lMAP,NTT,KfNT,KMN,MX 
OIMENSICN  AREA<2),  ARhC( 2 ) 

NPR«t> 

on  300  1*1,2 

c  ***** 

f  1-1  is  the  target 

C  -2  is  the  projectile 

c 

r  akea  is  the  initial  area  of  each  CELL  !N  square  inches 

C  ARHU  IS  THF  PROCUCT  AREA  *  DENSITY  IN  CLB-MICPnSPC-2/ IN-2) 

C  RHUT  IS  CNE-THIRD  RH0  IN  LB  MICROSEC  INCH  UNITS 

C  TWMU  IS  ThlCF  THF  RIOIOITY  MODULUS  AMU. 

C 

C  ***** 

AREA!  I  )  *  CRH)  *  ozm 
ARHOUI*  AREA!  1 )  *  Rul'd)  *  1.0012 
RHOTIII*  RH0m/3. 00-12 
TWMUIII*  2.0  *  AMUin 
?ro  CONTINUE 
c  ***** 
c 

C  NOW  oRID  IS  CREATED  FOR  BODY  l 
C 

C  ***** 

zcur  *  r.o 
jjj  *  jpin 

JJ2  »  JMI2I 
KKK  *  KM  III 
DO  7777  Ml, KKK 
1)0  7777  J*!,J4J 
7777  1PTAIK,J)«0 

00  326  Ml, KKK 
RCUR  *  C.r 
00  325  J«1,JJJ 
AR  ( K , J  i  *  RCUP 
AZ  ( K  ,  J  )  *  ZCUR 
ARH  ( X  ,  J  )  *  RCUR 
AZH  ( K , J )  *  ZCUR 
ARDH( K  ,  J )  «  ROCTU) 

AZOHIK, J)  «  ZDOTI 1 ) 

IF  1 .1-1)  3C4,3‘,4,302 

,  3 02  IF  (K-l)  204»3''4,3*3 

303  AA( K, J (*AREA( 1  I 
AV(K,J)*l.C 
AP I K, J  )*C  ,0 
ATRRI K, J)*0,0 
ATZZI K, J)*0,0 
ATRZIX, J  )»C.O 
ATTTI K,  J) *C«0 
ASRR (  K »  J ) *0.r 
ASZZCK,  J  1*0,0 
ASRZ(K,JI»0,'> 

ASTT (K,  J)*C«0 


46 


47 


L/  U  O  C  U 


n-5  ip  *  11 

AZDH(K,JI*2D!NT 
GO  TO  324 

320  IP  *  7 

GO  TO  324 

32 1  IF  (J-JJJJ  323*322,322 

322  IP  *  10 
GO  TO  324 

3? 3  IP  *  l 

324  IF  (IPTAIK,J».FC.ni  !PTAIK,JI«IP 
RCUR  *  RCUR  ♦  DP  (II 
3?4  CONTINUF 

WRITE  <NPR,70r*|  K,UPfA(K,KK  ) *■  KK*l« JJJI 
ZOJR  *  ZCUR  ♦  HMII 
■  32*  CONTINUE 

£«•«**« 

C  NOW  COMPUTE  CEIL  PSEUDO-PASS  ZM 
r ****** 

no  40 r  j*  ?,  u j 

TEMP*AP  C 1 , J) 

4'V  ZM(l,J>  *APH(im*(TrMP  -0.5*DRIliI 
***** 


NOW  GRID  FC«  BOOT  2  IS  CREATED 
***** 

ZCOR  *  ZCUP  -  DZ(I) 

KKK  *  KM (2  I 
JJJ  *  JM2) 

00  349  K*l , KKK 
RC'JR  *0.0 
00  348  J*1,JJJ 
8R  ( K  « J  )  *  RCUR 
DZ  <K,J>  =  ZCU» 

(>RH  ( K  ,  J  )  *  HflJR 
3ZH  (K,J»  *  ZCUP 
BRDHI K , J)  *  PrLT(2) 

RZl)M(K,J)  *  ZDOT  (  ?  ) 

IT  (J-U  33«,?V',3?p 

328  IF  (K-U  ?3'V33<',3?9 

329  i)A  ( K,  J  )  *  AREA!  2  ) 

(3V(K,  J)  *  I.-' 

8°  I K,  J  l*P  »r 
GTRRCK, J)*C.O 
‘3TZZ  (  K*  J  )*C.O 
3TRZIK, 

HTTTI*,J)*C,0 
8SRR  <  K » J  J *c.r 
OSZZIK.J  )*0.0 
RSRZIK.J)  *C.n 
8STTIK, J)*C.O 

r  ***** 

C  NCW  THF  PRIM  CONDITION  CnnE  IS  OF  TFRMINEO 

C  ***** 

ST'*  IF  ( J-l)  311  ,33!  ,334 
3>«  IF  (K-l)  332,38?,-»3-> 

»:<?  ip  =  13 

BZOHI K , J I *ZD I  NT 
GO  TO  34 T 

33 »  IF  (K-KKK)  334, <34, 334 


JH  IP  *  4 
00  10  347 
333  IP  *  ? 

GO  TO  347 

336  IF  (K-JI  337.337,340 
’37  IF  (J-J.)J)  330,’3ft,33« 

73P  IF  44'' ,339,6*0 

.3  34  IP*  10 

BZOHIK, JJ*/0INT 
GO  TO  347 

f)4P  IP»Q 

3ZUHCK,J>*?niNT 
GO  TO  347 
33^  IP*I 1 

HZOHI  K»  J) *201  NT 
GO  TU  347 

340  IF  ( K~KKK )  344,341,341 

341  IF  (J-JJJ)  343, 34?, 342 
34?  IP  *  4 

GO  TO  347 
34’  IP  *  7 

GO  TO  347 

344  IF  (J-JJJ)  346,345,345 
’45  IP  *10 
GO  TO  347 

346  IP  *1 

347  IPTBI K , J )*  IP 
RCUR  *  PCUR  *0R  C? I 

344  CONTINUE 

WRITF  (NPR.7000 )  K, ( I PTB(K,KK ) ,KK*I , JJJ ) 

2CJR  *  ZCUR  *C7(2) 

340  CONTINUE 
C  ****♦. 

C  NOW  COMPUTF  CELL  PSEUOO-MASS  7M 

C  ***** 

on  401  J-2.JJJ 
TEMP*BR  ( 1 » J  ) 

401  ZM( 2  » J  )  *  ARHli  ( 2  )  *  (TEMP  -',.F*0R(2M 

***** 

INITIALISE  CYCLE  COUNTF R ,  TIME,  TIME  STEPS, WAVF  BOUNDARIES 

***** 

ICYCL  *0 
TYME  *C 
OTN  *OTF 
OTHM  *OTH 
KMN  *K  INT  -3 
KMX  *  4 

JMX  *  JM2)  ♦  ? 

IF  (  JMX.CT.JMd)>  JMX*JM(  1 ) 

7000  FORMAT  (*C*,20I5> 

RETURN 

END 

SUBROUTINE  LETGC 

IMPLICIT  PEAL  *4  ( A-M.C-Z , I) ,  INTFGCR  (J-Nt 

COMMON  /CCM1/  AR  (24,46),  ARM  (24,46),  AR0H(?4,46), 

1  A2  (24,46),  A2H  (24,46),  AZDM(24,46), 

?  AA  ( ?4 » 46 )  ,  AV  (24,46),  AP  (24,46), 

3  4TRR ( 24  ,46 )  ,  AT2Z(?4,46>,  ATRZ124,46), 
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i 


ATTT<24,46) , 
ASPR<?4 ,46), 
ASTT(24 .461 


4SZ?  (24,46)-,  ASR7I24. 46). 


COMMON  /CCM2/ 


COMMON  2CCM3/ 


SP  (46. 
B?  (46. 
BA  (46. 
HTPR ( 46  * 
BTTT <  46  » 
RSRR (46 « 
HSTT(46  » 
TR  (24 


61  . 
61  • 
6)  « 
6)  t 

o!  . 
6  i  « 
6) 

.46)  . 


BRH  (46. 
8ZH  (46. 
9V  (46. 
0TZZ(46, 


6). 
6)  . 
6). 
6)  . 


8ROM! 46*  6). 
BZDH( 46.  6). 
BP (46.  6). 
BTRZ(4o.  6). 


BS7Z(46.  6).  BS*»Z(46.  6)  i 


TRH  (24.46).  TRnH<24,46). 

(74,46).  TZH  (24.46).  TZ0rt!24.4fc). 

SR(4fc),  SZ(46),  TSP.  ( 46) ,  TSZ(46),  ZM(?,46»,  SPT.^(S), 

AU2).  *2(2).  A)  (? ) .  A4(2),  A(FN(»),  AMAT(2>,  AM')(2).  *»••»■»)* 
BK( 2) ,  CAPE(  2)  *  COOT(?>,  f»R(2).  OZ(2).  RAO(?),  71 7U 

„  RHOG(2),  RHOT ( 2) .  TWMU(2),  YT(2),  ZOOT(2).  ZV(2).  VrON(2). 

7  CQ.DTN.PTH.OTHV.DTMf N»201NT»TVME»CH€Kn#FACTR 

COMMON  /  ICCM1/  |PTA( 24.46) 

COMMON  /ICCM2/  tPTB(46,  6) 

COMMON  /ICQW3/  ISHP(2),  JM(2).  KM(2).  KMIO. 

I  (RUN, 1CT1. (0T2. (DENT , (MAX, 1PRNT, ICYCL , IMAP.NTT.KINT.RMN.KMX 

JJJ«JM( 2 1 
00  350 

AR ( KM ( 0 . J )»BR ( 2 » J ) 

ARH  (KMIO»J)«BRH  (2.J) 

AR0H(KMIC.J)*BP0H(2.J ) 

AZ  (KMIC  » Jl *PZ  (2 . J  > 

AZri  (KMI0.4I-BZH  (?,J) 

AZOH(KP  !C.J)*BZ0H(2.J ) 

ATRR(KMID.J)*BTPR(2.J) 

ATZZ(KMIC.J)“BTZZ(?» J J 
ATRZ(KMIC.J)«BTPZ(?.J) 

ATTT(KM10,JI«BTTT(2,J) 

ASRR(KMIC,J»»»SPR(2,J> 

A5ZZ(KMIC,JI*BSZZ(?.J» 

ASRZ(KM10.J»»BSRZ(’.J) 

ASTTUMtC.J)«BSTT(2,J) 


AA 

AV 

350  AP 

**•* 


(KM!OsJ)»BA 

(KMIC.J)-BV 

(KMIC.J)-BP 


( 2 .  J  ) 
(  2 .  J  ) 
<2 .  J ) 


FINO  THE  SMFEP  BCONCARICS 


*4*4 

IF  (ZOCT(U)  300,301,300 
3C0  KMN  *  I 
JMXsJM(l) 

GO  TO  3f3 
IF  (KMN-1) 

KMN  «KPK-1 

if  uMx-jmn 

JMX  »  JMX  *1 
IF  (KMX-XM(2)) 

KMX  »  KMX  ♦! 


301 

302 

303 

304 
303 
3^6 

*44* 


3CJ.303, 1"? 


304,303,305 

306.307,307 


CHECK  IF  THE 
LIMIT  CF 


TIME  STEP  HAS 
10  PERCENT  OF 


ATROPHIED  BELOW  THE 
INITIAL  TIME  STEP 


ACCEPTABLE 
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nmno  o  o  o  r>  r>  o 


C  ♦*** 

!;7  IF  (DTHh-QTMlN)  10B,3C4,309 
3.C  1MAX  *  ICYCL  ♦  1 
GO  TO  3X0 

*«** 


N  N+l/2 

PINO  NEXT  DFIT  AND  OELT 

**** 

3^fl  OTN  «  (0TE+0THKI/2.0 
OTH  «  CTFW 
CHEK[)«0TH» 

*  +  ** 

INCREMENT  THE  CYCLE  C CUNT FA  AND  THf  TIME  VALUE 

«*** 

*»C  ICYCL  *  ICYCL  ♦! 

TYME  *  TYMF  ♦CTH 
RETURN 
E  NO 

5!lpMcIrLr«V,,If,nELZ’DeL2H'  xlst.kfnd.jh?) 

IMPLICIT  BEAL*fl  I A-H, C-Z , t) ,  INTFGFR  ff»wi 
COMMON  /CCP1/  A»  •-  *  INI  l-l,  FH  (I-N) 


CJNMUN  /CLF2/ 


4»  (24,46), 

A/  (24*46), 
AA  (24,46), 
ATp«r»4,A6)  f 
ATTT(24,46) , 
ASRR { 24  »46 ) , 
A  ST T(  24,46) 
OR  <46,  6)  , 
BF  <4  of  o )  t 
BA  <46,  6), 
0TRR ( 4o ,  6 ) , 
OTTT 1 46  ,  6 )  , 
BSPR(4b,  61  , 
HST  T(  46  ,  6) 

T°  (24,46), 
T7  (  ?4 ,46)  , 


ARH  (24,46), 
AIM  ( 24,46) , 
AV  (24  ,46 ) , 
ATZZ (24,46) , 


ARDH( 24,46 ) , 
A?HH( 24,46) , 
AP  (24,46), 
ATP Z I  24,46) , 


ASZZI24  ,46) ,  ASF 2 ( 24,46 ) , 

3RH  (46,  6) ,  HR DH( 46 ,  6), 

BZH  (46,  6),  BZOHI 46,  6), 

BV  (46,  6),  BP ( 46 ,  6), 

BTZZt46,  6),  BTRZI46,  6), 


HSZZf4o,  6)  ,  BS»r(46, 


'c™;"k;scvxTM^s;:!;;^r’TvME,c"Fm’F«r':  7,1  vrM'',>' 

CCMMON  /  JCCM2/  f PTB ( 46  ,  6) 

COMML  N  /ICCM3/  ISHP ( 2  ) ,  JM 12).  KM<2> 

lKl?K?:i?sU50T?Vf0,CNr,I'iAX,rPRNT*ICY^  ’  ^^.NT;,Kf.MT,KMN.K^X 

00  300  K«KST»T,KENC 
Du  .100  J*  1 , J(’2 
PZ  (K, J  )  =  EZ(R, J  I  ♦  OELZ 
xc r  BZH(K,J)  *  BZH( k , j )  ,  rei/H 
RETURN 
ENO 

*  SUHROUT  [NE  FMC4^’^S,JPT,KUX,F,Z,RH,ZH,RnH,ZOH,TRR,TZ7,TFZ,TTr, 

oMF,is!aRFAL*5'  n~*n 
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1  RUUX.i),  ZCKUX.l),  RU(KUX,l),  ZH(KUX,t>,  ROHCKUXtl  I  , 

2  ZOH(KLX.l),  TRRIKUX.i),  TZZIKUX.l),  TRZ(KUX.l),  TTmuX.il, 

3  SRR  (KUX. 1) ,  SZZIKUX.lt,  SRZUtJX.l),  STmUX.lt,  P(KUX,1> 

COMMON  /CCR3/  T R  (24,46),  TRH  (24.4ft),  TROH(24,46>, 

1  TZ  (24,46).  TZH  (94,4ft),  TZDH(?4,46), 

2  SR (46) ,  SZ (46) ,  TSR(4ft),  TSZ(46),  ZM(2»46).  SPTKf  5) , 

4  *1(2)  ,  >2(2).  A3(  2) ,  A4(2),  ALEN( 2) ,  AHATI 2) ,  AMU(2),  ANUl  ?)» 

ft  BK(2  I ,  CAPE(2I,  COOT (2),  ORI2),  OZ(2),  RAO(2),  RHO<2),  RDOT(2) , 
ft  RHUG12),  RHLT  (2) ,  TWKU(2),  YT<2),  20UT(2),  ZV(2)t  VCON<2>. 

7  CQ.DTN.CTF.OTrtW.OTMIN.ZQINT , TVMF .CHEKD.F ACTR 
COMMON  /ICCR3/  ISHPI2I,  JM(2I,  KH(2),  XMIO, 

1  (RUN  , 1CT1.IDT2. IOENT «!MAX, I PRNT, ICYCL  » IMAP, NTT »KINT »KMN»KMX 
NPR*6 
KCGL*KK 

IF  <KCCL<EC<1 )  RETURN 
Dll  302  J*2,JL$T 
SHR*TRZ(KK,J) 

IF  (OABS(SftR)-SHS)  3*2,300, 300 
yy*  ip  *  t pt( kk , j i 
JMl*J-l 
JPl*J»l 

IPL»IPT(KK,JN1) 

IF  ( (  (  IP.CT.13).AN0.UP.LT.20t>.0R.<(  IPL.GT.13)  .AND.  II  PL.  LT. 20) )  1 
1  r,0  TO  302 

WRITE  I NPR , 1C 0 )  KK,J 

100  FORMAT  ( 1F0,  •*<■**  SHF  AR  FRACTURE  ****  AT  K«*  ,  lft,‘.X,  •  J* •  , Ift /) 


TR  (KCLL.J 

) 

s 

R  (KK, J 

TZ  (KCCL.J 

) 

* 

Z  ( KK. J 

TRH  (KCrt,J 

) 

* 

RH  (KK, J 

T/H  (KCCL.J 

) 

* 

ZH  (KK.J 

TROrl(  KCLL.J 

) 

*» 

ROH( K  K, J 

T/OH( KCCL  , J 

) 

3 

ZOHIKK.J 

IF  ( IP.EC.R)  GO  TO  ') 

IF  (lP.EC.l)  GO  TO  l 

if  iip.fc.7)  r.c  tc  y 

GO  TO  3*2 
I PT  (  KK , J )  *  14 
IPT(KK,JP1)*21 
GJ  TO  330 
l  IPT(KK,J)*lft 
IPTIKK.JPl)*?? 

GO  TO  33C 
y  IPT(KK,J)*lft 
!PT (KK , JF 1 1  =  ?  ? 

33"  CONTINUE 

DO  30rc  11*1,? 

JJ*J*U-l 

KM1 *KK- 1 

KMINl=KPl 

JM1*JJ-1 

TRRW*  TRR  I  KK  ,  J.)) 

TZZW*TZZ(KK,JJ  1 
TR  Z  W*  TR  Z  (  KK  ,  J  J) 

TT  TW*T  T  T (KK,  J J  ) 

IP*IPT (KK.JJ) 

Rl*TO IKK,  J J ) 

Z!*TZ ( KK,  JJ ) 

( PLEF  *  I PT (KK-l  ,  J J I 

(F  ((lPLFF.GT.13).ANf).(lPLFe.LT.?Al)  :)0  TO  51*’ 
GO  TO  51C2 
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51*1  Z2*TZ(KK-t»JJ) 

R2*  TR ( KK~ l ,  J  J ) 

Gu  TU  51C3 
510?  Z2«Z(KK-1,JJ) 

R2«R(KK-l,JJ) 

5103  CONTINUE 

IPR=IPTIKM1,JM1) 

IF  C(IPR.GT.!3J.ANn.< IPR.LT.20) S  GU  TU  5019 
GO  TO  5018 

50I<*  Z3  *  TZIKPINl, Jvl  1 
R3  »  TR<KFINI, J«l> 

GU  TU  5«?/' 

5018  Z3=Z ( KM1 ,  JM1  ) 

R3*R(KMI,JM1 » 

5020  IPR  =  lPT(KK,JMl  ) 

IF  ((  IPR.GT.m.ANO.  (  IOR.LT.?^))  GO  TO  5021 
GU  TO  502? 

5021  Z4*TZ(KK,JJ-l I 
R4  =  TR  <  K  X  »  JJ- 1 ) 

50  TU  5C23 

5 32?  Z4=Z<KK,JPl> 

R4*R(KK, JM1 ) 

50?3  CONTINUE 

CALL  FRFESl  ( 1 ,R1 ,Z1 ,o?,Z2,R3,Z3,R4,Z4,TRRW,TZZw,TRZW,TTTw,SRPW, 
lbZZW,$0Z).,STTW,KK,  !P,  Pw) 

TRRIKK,  JJ)*TRHW 
T7Z(KK,JJ)sT/Z2 
TRZUK,  JJ)=TRZW 
TTTIKK  ,  JJ)*TTTW 
P ( KK , J J ) *PR 
SRH(K5,JJ)=SRRk 
SZZIKK, JJ)*SZZa 
SRZIKK, JJ)=SPZW 
STT(KK,JJ)*STTrf 
30('**  CONTINUE 
302  CONTINUf 
RF  TURN 
FN') 

S'JRRO'JTIAE  FRACTN  (  TRP  »TZ  Z »  TP Z  »TTT  ,K  ,KMX  ,R,  Zf  RH,  ZH, R  OH,  ZOH,  J 1 , 
ITSTR,SRR  ,5ZZ,$R7,STT,P,IPT> 

IMPLICIT  REAL  *8  (A-H,0-Z,t),  I NTFGFR  (l-N) 

OIMENSICN  I ° T (  Ktjx *  1 1 » 

1  R(KUX,1),  ?(KUX,tl,  RhlKHX,!),  ZH(K')X,1),  RI)M(K'JX,1  I  . 

<:  ZDHI  ALX  *  1 )  ♦  TRR(KUX,1),  T7  Z I  KUX,  1  ) »  TRZIMIX,  » ),  TTT(KUX,1), 

3  SRR  (KL'Xil  )  f  SZZ ( KUX, 1 )  *  SPZ(KUX,1),  STTIKUX.I),  P(KUX,n 

COMMON  /CCM3/  T  R  (24,461,  TRH  (24,46),  TR0H(?A,46>, 

l  r 7  (24  ,46),  T  ZH  174,46),  T7 OH! ?4, 46 ) , 

?  SR ( 46  )  ,  SZ ( 46 ) ,  TSR(46),  TSZ ( 46)  ,  ZM(2,46),  SPTM(5), 

4  A 1  (?)  ,  62(2)  .  AM?),  A4(7),  ALFN(?>,  AMAT12),  AMU!?),  AN'H?), 

5  BK  (  2  )  ,  CAPF  (  ?  )  ,  f.nOT(?),  ns(’|,  nZ<2),  PAr>(?),  RHl)  ( ?  ) ,  ROUT!?), 

6  RH0G(2),  PHOT  (  2  ) ,  T  W  M»  J  ( 2  )  •  YT(2),  Zn0T(2),  Z'/(?l,  VCONCM, 

7  CQ,r)TN,rTH,rtTH«,OTMI  6,201  NT,  TVMr,CMfKD,FACTR 
NPR=6 

KPl=Kd 

00  3">(r  J=2, J) 

STR=  TRP ( K  ,  J ) 

STZ=TZZ(K,J) 

$  TR  Z*  TP  7  (  K  ,  J  ) 

FCUN=0.5* ( STR+STZ ) 

T0UN=0.5* (STZ-STP ) 
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IF  (OABSITCONM  430,430,420 
4? 9  SBC*$TRZ/TCCN 
Prtl*OATAN(SBCI 
GO  TO  431 
439  PHI*O.C 
431  ALPHA*0SlMPhI) 

BETA=OCCS  (PHI ) 

TZZP*ECCMTCON*HETA  ♦STRZ*ALPHA 
TRRP=-TZZP*$TZ*STR 
PKI*90,0»PHI/3. 141S928 
IF  (rZ7P«TRRP>  3C2, 393,303 
30?  IF(TPRP)  312,312,314 
314  tns*tppp 
GO  TO  322 

3“M  IF  (TZZP)  312 ,312,374 

3?4  tns=tzzp 
GU  TO  322 

32?  !«=  (  TNS-TSTRI  312,326,326 
??6  JM1= 1-1 

3?0  WRITE  ( NPR » 10* )  K, J , PH l 

IC9  FORMAT!  IH9,‘*«*«SPAU**F*AT  <  =  •  ,  1 5 , 5X » • J* • 1 5 , • ANGLF  = ' FR .?  ) 
IPT(K»J)=1° 

TRR!K,J»»C.9 
TZZ<K,J  1*0.9 
TP.ZIK,  J)«0.0 
PW=-TTT(X,J)/3.3 
POC,J»  =  PU 
SRRIK, J»*PW 
STT<K,J)*-?.P*PW 
SRZ ( K, J  1*0.0 
AZ  Z  I  K  *  J )  *Pl» 

TRR(KP1,J)=0.0 
TZZIKPl  ,J1*3.9 
TKZIKP1,  J)*P.O 
PW=-TTT<K,J)/3.0 
PtKPl, J)*PW 
ARR (KPl » J  J*PW 
AZZIKPl  ,J)=PVs 

ARZ(KP1  ,  J)*9.o 

STT(KP1,J)=-2.0*PW 

ZO=Z I K , JP1 ) 

RO  =  R  (  KPl  ,  J) 

ZA*Z(K-l,JPl) 

RA=R(K-1, JP1 I 
ZB*Z<  K-l  ,J) 

R8*R I K- 1 » J ) 

ZC=Z( K.J) 

RC*R(K, J  I 
IP  =22 

TRP  R*  TPP ( X , JP 1 ) 

TZZW*TZZ(K,JPU 

TRZW=TRZ(K,JP1) 

TTTr*1  niKfJPU 

CALI  FRFESO  (l,RO,ZO,RA,ZAfRH,ZB,PC,ZC,TPRW,TZ/w,THZW,TTTW, 
l  SRRW,  SZZV>,A»Zl»,STTW,K  ,I<»,PW) 

TRfUK,  JP1)=TRRW 
TZZIK,  JP1)=TZZW 
rRZIK, JPl)=TP/W 


I 
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r»  o  *■> 


tttu.jpi  j«tttw 

P<K,JPU=P* 

SRR{K,JP1)«SPRW 
SZZ(K, JP1 )*SZZW 
SRZ(K» JP1  l»SR2i( 

STT (K , JP1 )»STTW 

IPT(K,J*1)»18 

JPi=J*I 

IPT (K» JPl)=22 

IPl*IPT(KPl,J) 

If  (IPI.EQ.22)  IPT ( KP  1, J ) =20 

TRIKt JJ=20 

TZ(K,J)«Z(K,J» 

TRH(K, J)«RH(K,  J) 

TZH(K, J)«ZH(K,J) 

TROH(K,J)=ROH(K,JI 

TZDH(K,J)*ZOH(K,J) 

TR ( K  »  JM ) )=R(KfJM]) 

TZ(K,JM1)  =  Z(K,  JMl  ) 

TRHIKi  JPII*RH(K,JMH 
TZH(K, JM1 >*ZH(K,JM1 ) 

TROKU,  JP1»=RCH(K,JM1 ) 
rZDH(K,JMl)=ZDH(K,JMl > 
il  2  CONTINUE 
TCC'1  CJNTINUE 
RETURN 
END 

SUBROUTINE  GUTT  ( IR  1 1 PT  ,R  ,  Z ,RH ,ZH tRDH , ZDH , TRR f TZ Z » TRZ i T TT , SRR i 
!  SZ7, SRZ,STT,A,V,P,KUX,KMIN,KMAX,  UMAX) 

IMPLICIT  RFAl*fl  (A-H,0-Z,t),  INTEGER  U -N ) 

COMMON  /CCM?/  TR  (24,461,  TRH  (24,46),  TRDH(24,46), 

1  TZ  ( 24  f  4h | f  TZh  (24,46),  TZDH(24,46), 

2  SR (46 ) ,  SZ (46 ),  TSR ( 46 ) ,  TSZ(46),  ZM(2,46),  SdTM(5>, 

4  Al  (  2 ) ,  *2(2) *  A3 ( 2 ) ,  A4  (  2 ) ,  ALEN(2),  AMAT ( 2 ) ,  AMU(2),  AN0I2), 

6  9K(2),  CAPE  ( ■*  I ,  COOT!?),  0R(2),  0Z(2),  RAD(?(,  RH0I2),  RD0T(2), 

6  RHOG  <  2  ) ,  RHr ''  ( 2  ) ,  TWMU<2),  YT  (  2  I  ,  Z0PTI2),  7V(2>,  VrON(2), 

7  CU,DTN,CTf,DTHH,nTMIN,ZOINT,TVMP,CHFKD,FACTW 

CUK-tUN  / ICCM3/  ISHP(2I,  JM(2),  KM(2)f  KM [ D , 

1  I  RUN, I  DTI, ID T2, IDFNT  ,  l M AX , I PRNT , ICYCt  ,  (MAP , NTT ,K INT ,KMN, KMX 
D1 AFNSICN 
l  I  PT  (Kl)X,I )  , 

1  R ( KUX,  l  )  ,  RH ( KUX  ,  1  )  ,  RDH(KUX.l), 

2  Z  (KUX ,  l  ) ,  ZHIKUX,  I),  ZnH(KIIX,l), 

^  TRR ( KUX ,1 ) ,  T / 2  < KUX ,1 ) ,  TRZ(KUX,l),  TTT(KOX,l), 

?  SRR  (  KUX,  1 ) ,  SZZIKUX, 1 )  ,  SRZ  (  KUX  ,  l  )  ,  STTUllX.l), 

5  A ( KUX , 1 ) ,  V ( KUX , 1 ) ,  P ( KUX , ! ) 

NPR  =  6 

SHS  =  AI (  IR) 

**  ** 

INITIALISE  TC  ZERC  VARIOUS  VARIABLES 

C  **** 

TRR A  =C.O 

TRRB  =C.C 

TRRC=O.P 

TRRO=O.C 

TZZA  =  O.T- 

TZZB=0.r 

TZZC=0.0 


55 


tr2a*o.o 

TR?H=0.a 

TRZC**''./' 

T82D=0.C 
AC  ON A  =  C . a 

AC JNS=A  .A 
ACONC=C.C 
ACONO*C.r 
i)CoNA  =  0.0 

HCon8*a.a 

lie  ONC 

arcwso.c 

PCuHA*C.C 
PCuNH=f  .' 

PCONC =C.r 
PC(JNO  =  r  .A 
!J2=0.0 

n  3*0. c 

234=0.*' 

24 1  *0 . 0 
PU-sf'.r 
R?3=a.a 

R34=0.C 
4  1  *  n  .  ** 

Af)  1001  K=  K**  I  A  »  KMAX 
Kill  .41  =  K-l 
KPLSl  =  K  ♦  ' 

KMl  =  KMM 
KC  HR*  * 

API  =  KPLSl 

4*"  l  00  3  3C  A  J  =  1  ,  Jf* AX 

JMl  =  J-l 
.IPl  =  )♦! 

IP  *  IPTIK.J) 

I P  (K.CT.l)  IPLPT.IPT  (KMIMU  JJ 
Si  3  1  t  T  A  =0 
IT3=C 
ITC=0 

ITO  =  n 
11$-' 
i  r  iA*o 
IT'i.js' 

|TOC»() 

8'"i  r,j  Tu  (SfCl  .si  oi  ,5Pnl  .SAOl  ,  ,80^1, 5-04,5r"4,*,A0l  ,50^1,51' Al, 

1  *>r01.  ‘jA".  ,  5' A  l  .Wl,  S '('4,5^34,5'''' l  ,SrA4,  SAP  1  ,S004,  KOOl  . 

'  SO  34, 5004,4  00,  *0*4  ,kao4,Kaa.2»SC'A2).  IP 
'I  Rl  -  R(KPl.J) 

21  =  2IKPI.J) 
r.u  to  scr 4 
SC  IPKMPTIKPl,  J  J 

IT  (UPR. GT. 13). ANO. (  IPP.LT. 2C.M  GO  TO  SA<m 

go  rn  sfpt 

SA'J  »l  =  TRIKPISl.Jl 
21  =  T2IKPIS1 ,J) 

■jaA4  Cl 'NT  IN' IF 

)  TO  (5005,500s, 5005, KOOK, 5005, SA'^.AOASjSCAA.Kf  ak  ,  5.}''At  S'  'S, 
l  SCCS.5AA*  ,  5^ A  »5A.1A*4A0(,  f  «;aos  f  ^nr>is,5nrt,  , 

i  5015, *005, 5005,5005,5006, 5224, 5226) ,  IP 

sc  's  it?  =  m kcijr  »jpi  ) 
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Z?  *  Z<KCUR,JPl) 

522E  IPR*IPTIK»JP1  } 

IF  <<IPR.GT.13).AND.( IPR.LT.2n))  GO  TO  5227 
GU  TO  50C5 
5? 2 7  R2»TR(K,JPl> 

Z2«TZIK,JP!) 

5f'6  CONTINUE 

GO  TO  (5fC7#5Pr7,50n«5rt07,5''l''»5^tC'f5'jn7,5rtr'7,5O07t5<'f'7,5rnr, 

1  5007 ,5007,5007, ‘5007, 500^, 5010, 5007, 5007, 5000, 5007,5007, «007, 

2  5''3fl,50C9,5CCft, 501 0, 5010, 5009  ),tP 

5f'7  R3*R (KM1 , J ) 

Z3»Z(K*l,J) 

Oh)  TO  (5036, 5*34)  ,18 

5036  IF  UK.EC.KMINl.ANO.lK.GT.in  GO  TO  5036 
5"34  Co.*TU(jE 

GO  TO  5*10 
5016  i«(J)sR3 
SZ( J)*/3 
GO  TO  501'' 

6°  "'<?  «3=TR(KP1,J) 

Z3*TZ(KP1,J) 

Gil  TO  5010 
•iC'H  IPR  *1 P  T ( KM1 »  J ) 

IF  (UPR.GT.131  ,ANO.(  IPR.LT.20>)  00  TO  60*0 
GM*  R3*RTfivo 
Z3=ZTtPP 
6 *  1  "  CONTI NUF 

GO  T(1  (  5011, 50 1 A, 5016, 5*14,5*11 ,5011,5011,501  1 ,5911 ,5*11 ,6f  11 , 

1  5011, 5016, 5011, 50 11, 5011, 50 14,5011, 6012, 5013, 5013, 6013, 

2  601 3, 5*14,5" 14, 501 4, 5*14, 5014, 5*  13),  IP 
611'  1  °Rs I P  T ( K , JM 1 ) 

IP  ({  IPR.GT.l  3).  AMO.  (IP'S.  IT.  20»)  Go  T0  I  B 
6-12  R4  *  R(KCLR,J31I 
Z4  =  Z<KCtlR,jMll 
Go  TO  5' 14 
VH’  K4  =  TR  { K  ,  JM  1  ) 

r*  =  tz<k,j»i  j 
6-14  continue 

GJ  TO  (4COO,4001 I , 1R 
4C"0  IF  u-ll  4CC7, 4103,4071 
4 IF  (K-KMIN)  4'’',,,477?t4'li'» 

40t),»  SRI  .1)  *B  I  K*1  ,  J ) 

szi  j)»;ikpi,.j) 

4^  ■>3  r  ii^TiNiir. 

7-6  GO  TO  (  1,2, 3, 4, 5, 6,7,8,9,r,U, 12, 13, 14, 15, 16, 17, 16, 19, ?i, 21, 22, 

1  27, 24, 25, ’6, 77, 78, 78}, jt» 

1  I T  A  =  l 
ITB=l 
ITC-l 
I  T0»1 
GO  TO  7 C  4 

2  I  T  5  *1 
I  TA  *  1 
I  TH*  1 

G!)  TO  70? 

1  I T  S*  l 
I  T  A*  1 

Go  TO  7  C 1 
4  I T  S=  1 
I  T  3s  l 
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ITHA«1 
GU  TlJ  702 
S  I TA*l 
I T  0*  1 
l  T  l)C  *  l 
GU  TO  704 

h  iro»i 

iruc»2 

uO  TO  70  4 
7  !T0«1 
!Tf.*l 
IT(1A*1 

go  rn  701 
a  ITC=1 
I  TOtl*  l 
Gu  TU  7C3 
o  ITB-l 
I  TC  =  l 
I  TO*  l 
IT  0  A*  1 
GU  Tu  704 
>  1  I  T  C  *  1 
!  r  o*i 
l TUA* 1 
oil  TlJ  704 
U  CONTINUE 
.0  TO  1 
12  CONTINUE 
Gu  Til  1 
'  ’  I T A* 1 
I  T  A*l 
I  T  S*  I 
GJ  TU  7r2 

14  O.oNT  I MJO 
GO  Til  10 

15  CCNI I MJF 
Oil  1U  l** 

16  CONTINUE 
Gil  Til  8 

’  t  CUNT  I  Nil  C 
GO  TU  3 
1J  C.UNT  INI.F 
«JU  III  « 

1=  CUnT  I  M.F 
GU  TU  P 
20  CONTINUE 
Gu  TO  l 
7 '  CJNTINGT 
GU  TO  7 
2 2  CONTINUE 
GU  10  l 
?1  CONTINUE 
GU  TO  7 

<;4  CONTINUE 

I  TH  =  l 
ITOA*  ! 


GO  TO 

70.  ? 

2S  I T A*  1 

l  T  A  *1 

Gu  T  u 

7  62 

26  CONTINUE 
GO  TO  24 
?7  CONTINUE 
GO  TO  3000 
2*  ITA«1 
GO  TO  70! 
2$  ITA»1 
IT8*1 
ITIJ«l 
IT0C«3 

™4  Z*l  «Z4-Z1 
R41  *R4-E1 
TRRO  «  TRRA 
TZZO  »  T?ZA 
TRZI)  *  TPZA 


ACO  TD*  ACTIN'  6 
BC.UnD*  HCCM 

IE  (ITCC.E2.TI  G(1  T0  7?? 

*  I TOf-21  306,313,706 
-1'  IE  (ITCC-l)  763,701,703 
Tr-T  Z34.Z3-Z* 

R34*R3-R* 

TRRC  *  TRRO 
TZZC  *  TZZP 
TRZC  «  TRZB 
PCONC*  PCCISB 
ACOVC*  ACCNtt 
rtCvi.3C«  PCCNB 
ip  (irrp-i)  3^7, 3(3,707 
307  IF  (IT0H-2J  7C?,7d,707 
7*2  Z2J  «  Z2-Z7 
R? 3  *  R?-R3 


* r'  '  *  -  IPKIKIUP  ,JP1  ) 

TZZd  «  TZZ(KCU<5  ,JP1) 

TKZR  *  TPZ(KCU«,JP! » 

AN  *  A  ( KOUR  ,  JP1 ) 

AMASS  *  AN/Z*(  IR, JPl) 

PCONB  «  RHL(fR)*A:j/  V}KC'IR,JP1) 

AC GNU  «  TRZB+AMASS 

TTT(KCHR,JPln*AMASS 
’  *r  UTCA-IJ  76l,3l3,7Al 
7*1  Z 12*  Zl-Z? 

R12*  K)-R? 

TRRA«  Tor(kpi,jp(  ) 

TZZA*  T7/(KPl, J0(  j 
TRZA*  T»Z(KP( , jo( ) 

AN  *  A  (KPl.JPl) 

GC  TO  ( 306,311), I « 

3^  IE  (K-KIMI  71)  ,JK,3U 
■  1 n  AMASS  *  AN/ZM) 2, JP( | 

PC ON A  *  RHC (2  I  * AN/  V(KP1,JP1> 

GO  Tfl  312 

711  AMASS  *  AN/7M( (P, jpl | 

EC  OVA  *  RH)(  I<3  )  *AN/v(  KP1  ,  jp  1  ) 

312  AC'JNA  *  TRZA* AMASS 

OC  J.mA  =URRa-ttt(kpi,  jpi  ,  t,AtMSS 

31’  IF  (IP.GT.23)  GO  Tn  ^IP 
** 


o  n  n  n  ^  O  O  o 


EQOAT  tCNS  OF  MOTION 


**** 


f  ZN*ZCKCUR,J1 
RN«  R|KCLR,J> 

RDHF*ROHf  KCUft  «  J I 
ZDHF.znHIKC'JR.Jl 
GJ  TO  SC  17 
6  ZN*TZ<K,J1 
RN«TR(K„»} 

RDHF*TRCH{K,J) 

?f)HF*TZCHK»JI 
5C17  CONTINUE 

Prti*  UTA*PCONA  ♦  !Ta*PCflNR  ♦  !TC*PCQNC  ♦  lT0*l»CJW*»*l.',ni>/4.P 
RZOCN  *  fTK/<?*BMl> 

If  UTS-11  314,314,315 
314  RNR*r>.0 
RNH*?.'' 

KONit*-*.1' 

ALFA*0.0 
GO  Tti  318 

MS  TC-IH*  l.C/UTAMTBUTCMTO) 

ALFA*  TCCN*(ItA*AC^NA  ♦  !T3*ACUN8  !TC*ACONC  ♦  ITn*AC0N01 
BETA*  TCCN*(  I’A«8CONA  ♦  ITB*8CON«  ♦  ITOBCONC  ♦  ITI'-BCON'*  1 
RONH*  RPFF  ♦  F7CCN«UTA*ITRRA*Z12-TR7A*Rl?l  «•  ITR*(  TBPB*723 

1  -TI*7B*P2J1  ♦  |TC*ITRRC*m»rR2C*R34»  ♦  !TD*tTRR0*74l-T»Z0* 

?  R4111  ♦  l)TN«»ETA 

ROEF  *  RCM-*RTh 
RNW  *  RNAROEF 

IF  ( OABSIROFF  )-  l  ,00-07)  316,314,317 
716  RVV  *  9h 
*JT  Kf,M  *  mv,-»RNJ/V’ 

Ji»  ZUNH*  ZOFF  -P/OCN* U  TA*I TZZA*R12-TRZ  A*Z 1? 1  ♦  ITB*( TZZB*P’3 

l  -T*7B*Z73)  ♦  !Tr.*n7ZC*R34-TPZC*Z34I  ♦  IT0*(TZZt>*R4l-TRZ0* 

?  ZAlll  ♦  CTN*AUFA 

ZOEF  «zrAH*OTH 
INH  ■  ZN^ZCEF 

IF  <»ABSIZCFF)-l,00-07)  311,319,320 
31 9  ZNW  •  Zh 
3’"  ZNrt  *  <ZK  ♦ZNWI/?.'' 

IF  (ITC1  322,321,327 
371  l»H  1  *  * 

Gil  TO  «70 


**«* 


NOW  IMF  STRESSES  ARE  CCHP'ITCC 


»«  ** 


3,  > 


SI*  I 


*r  ?. 


zi*  ZNW 


Rl*  R.\» 

IF  ((  lRLFF.GT.tU.  ANO.IIPLEF.it.  ?r>)  GO  TO  81*1 

Gfl  TO  MT2 

Z?*TSZ(.n 

R2  *  TSR  ( J  1 

RH2*TRM!KI*1  ,  J) 


Zrt?*TZMK*l,J) 
R0H2  *  T  RHFt  ( K“1  ,J) 
ZOH2«TTrMK*l,  J> 


NOT  REPRO°UC,BU 


GO  TU  Flo  3 
Z2*SZ<  .1) 
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R?«SR<Jt 
RH2*RHOO*l,JJ 

2h2*2H( kp 1 1 J  > 

ROH?*ROMfKMt,J» 

ZDH2*ZnH<RM|»J) 

CONTINUE 
t  PR*| PT (KP1 • JMl  ) 

IF  (UPR.GT.t  JJ.AMl.I  IPR,LT.20|»  bl  TO  50l«J 
GO  TO  501S 

5015  23  *  TZIKPIM.JPII 
R3  »  TRIKMINI, JPt) 

ZHU  *  ZAH-TZHIKMINI « JM1 ) 

R«H  *  RRH-TRHI  «M!fl  I  *  JMl ) 

Zi)Hl3  *  ZONH  -TTOHIKMINl*  JMI  > 

ROM13  *  POKH  -T»DH<KMINt.,Mi) 

Gi)  TO  50  2" 

5015  ZWUPl.JMl) 

R3*RIKMl,JMt) 

Zul  3  *ZNF  -  ZH(KM1,JM1) 

RHl  )  »P.\|-  -  PM(K*1,JM1) 

ZDH13*/f5KH  -ZI)H(«P1» jmk 
RDHU*Pr>5E  -HCHtMlf J»l) 

5 o?o  Ipk=IPT(X,JM!) 

IF  ((I»«.  CT.H>.  AND.  (  Ii*».(  *.2'))  GO  TO  5021 
GO  TO  5'V* 

'VI  Z4«TSZ<JM! 

R4*TSR<J*1I 

ZH4?  *  TZHOVJMl)  -ZH2 
RH4?  *  TRHK,  JM1I  -RH? 

R0H42  =  TPrHlK.JMl)  -  RUH2 
Z0H42  *  TZriH(K,J"l)  -  ZOH? 

GJ  TO  50?  3 
5  V?  Z4  =  SZ< JP’  ) 

R4  *  SP(J“1) 

«U42  I®H(KC0W«JV1 J-KH2 
ZH42  •ZHIKClJDt JMD-ZM? 

Z004?*ZCM5CH«,  IMI)  -  ZflH? 

RD042  *ROH  {KC(J°  »  f <5 1 )  -  oom? 

50?  3  CONTINUE 
r  ***• 
r 

r.  COMPUTE  THE  APF4  AN"  KELA  T1 VF  VOl’.IMF  \r  \fW  ![wc 
C 

C  *«** 

AN  =  A  ( KCI.R •  J  1 
VN*  V(KCtH,J) 

Z42=Z4-Z< 

R42*R4-R? 

A 1 240.  f *(/**( Pl-t 2) -Z1*R42  0?*<R4-Rl  I  I 
A234=0.5*<Z4*fK?-O>«-Z?*<P)-P‘.MV3*R4?> 

AW=A124+A234 
AU=n, 5*  <  AW+AN) 

VW  =  »H|)T  (  I  R  )  *(A1?4MX4+F  UR2)*A2)4*<M3H»4«-R2)  )  /  LM  [rt.JI 
VH  =  ,>.  5->  ( V»  +  VN) 

DEL  VVS  ( VX- VN >  /  >/H 

IF  (jABSirFLVVI-1. 00-07)  >?).  )?3,  324 

’?)  DFlWO.f 

C  **** 

0 

C  CoH’OTF  THF  STOlb  lnC.RCrtFNTS 
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c 

r.  »* «* 

324  ECtJN  »0.«*DTH/AH 

EZZH  «  ECCN  *  CZ0H42«RHIJ  -  RH42*Z0H13I 
FRRH  «-ECCN  *  |R0H4?*ZH13  -  ZH42*Rf)Ht3l 

EKZH  *  ECCN  *<K0rt4?*RHn-RH42*RnH13-Z0H4?«ZHl3*Z0Hl3*ZH42> 

CTTH  «  HELVV-6 Z7N-ERRH 
SCONt  •  OEtVV/».rt 
SZZN  *  SZZIXC'JR.JI 
SRRN  *  SRMKCU*»J) 

SKZN  »  SRZIKC'IRtJI 

41 E A  «  C.M*EC0M<RI)H42*RHl3~RH42*RfJH13*ZnH42*ZH13-ZH42*ZnHt3) 
HE TA  *  2.0*SRZN*4t FA 

r.  *«■»« 

c 

r.  COMPUTE  NEW  STRESS  nCVIATCRS 

r, 

C  **** 

SZZW  «  SZZN  ♦  TWMOI  !  fi 1*1 EZ  ZH-SCCUf. I  *Rf  T A 
$RR M  «  SPUN  ♦  TWMU(IHI*<FRRH-SfMNC»-RCTA 
SRZW  *  SP2N  AM*'(  ? (3 )  *ERZM  ♦  (  S&RN-SZ  ZM I  *ALFA 
STTW  *  S*T  ( KCIJR  »  J  )  ♦  TWMUII  B)*(  ETTM-SCONC  J 

C.  **  ♦  * 

C 

C  CJMPUIF  THE  ARTIFICIAL  VISCOSITY  TERMS 

C 

C  *♦*» 

IF  (DELVVJ  32A.32S.325 
3? *»  OQNH  *  H.c 
GO  TO  3?7 

32‘  V'ljV  «  DEL VV/OTrt 

CORil  *  (  AH*P MO (  l  B  I/Vh  I  *(  CQ*VOPV  1**2 

r  *»  ** 

c 

C  USE  TmC  EQUATION  cr  STAT  =  TO  COMPUTE  THE  PR^SS'lRc  INCREMENT 
C 

r  *»•  •* 

•»? ?  na1*  *  -RKIIRMOFLVV 
pw  *  pikcur.jmoflp 

c  **** 

C.  TOTAL  STRFSSES-  ARE  CCMP'JTfD  NFXT 

c  **** 

TZZH  *  S7ZV-PW-CCNM 
TPRR  «  SPRM-Pr-CONH 
TRZW  *  SRZW 
TTTR  *  ST  Th-PV-l'CNH 

C  «♦  •• 

r. 

r  CALL  THE  FRFF  SURFACE  RC'JT  Irte  to  ADJUST  the  STRESS  vaiufs 
f 

r  t< 

CALI  FHMSCIlHfRltZl  ,R?  .  Z?,  »  3  ,1 3  ,R4,  Z4,  TPR’Wt  TZ  ZW*  TRZ  W.  T  TTW, 

1  SKR*tSZZW.SPZW»STTH,KtIPtPH» 
f  *r  ** 

C 

C  COMPOTE  IFF  PRINCIPAL  STRESSES 
r 

f  fr*** 

■>2R  KvJN  =  r.S*(T7/rf*TPpR> 

TCCN  *  r.  S*  (  T  Z  /Vt-TPRW  I 
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IF  IDAeSITCCNM 
J2>?  PHl=OAT*MT»ZW/TCfNJ 
GJ  Tu  33 1 
J3P  PH  1  #*' 

33!  ALFA  «Q$!MPHI» 

BETA  »CCrS(PFll 

U/l>  *  ECCN  ♦  TrON*BFTA  ♦  fR  7W*Al  FA 
tprp  *  -tzzp  *r;n,  *  trrw 

i»Hl  *<*'*.  P*PHl/?.l4li<;?B 

c  *<*♦* 

r, 

C  CHECK  FTP  TIME  STFP  STABILITY 

r 

C  *"  ** 

33?  RP  I  3*  1R1-R3J**?  ♦l/t-Z3l*’»? 

RR?4*  F42* R4?  ♦  Z4?*T4? 

IF  |RR!3  -  PR  24 }  333(333,334 

333  OFI.R  *I)SC°T(P:3I  31 
Gu  Tu  3  35 

334  itFLP  *OSC«T(P»  ’31 

3?.4  CHEKO«(nfLP/f.t:t:if <m  > *factp 
p»<3  CJNTIwor 

Oil  T J  (?3«,3*<M,  !M 
33ft  IF  <K-KMtM3  37,337,33q 
in  if  u-ima* 3ij,334 
HP  KT-K 
JT*l 

DTrtWxChFKC 

3<C  IF  <CnF*r-OTHt,l  34^  ,34'', 34! 

>4>*  KT*< 

JT  *  J 
IHT=I“ 

'Hmw  =  fhekt 

J4’  I f  HJCC  If.VfL,  IPRNTI )  351,342,351 

34'  IC  (  j- 1  )  34^,343,3«;r, 

3',i  if  ,  J  f»-  I  J  34ft  ,i-,4, 34ft 
<44  ||  (  K  -<ft  |  M  J4  ),  345,340 

r  ** *■  * 
r 

f  H'4.1  ll'F  MV  ftJC(ST  PCS!  T  JiJ  , 

c 

r  ** 

345  Mf-  -  TYPOrr.JT!  !  ) 

KRF  siCIMUr/'lM  1  |  +  '.S) 

(.XITF  ( NPR,  100  )  irYC.I.,TVMC,wr,KPF,')TH,rtTHW,inT,KT,JT 
34ft  ic  (  H-2)  .340,  347,349 
'4’  IF  IK-?)  144, >43,140 

34  p  W-U  TP  (f.Pl.t  31  |  <(NT 
34  f  XR1TF  <  4  PP  ,  1  ^  J  IR,K 
17 'P  -  0.0 
TPFP  i  ft.  ft 
PHI  =  *.r 
V  I*  *  ~.r 
TPR V  =  ft. ft 
TRZw  =  0.0 
Mrw  =  r.ft 

f  »«.'!* 

C 

ft  WP  l  T‘:  l  c  T  l'*r  (FS'll  TS 


0^C 
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n.Tnnnn 


C 


**#• 

3 SO  WRITE  INPR.IP3)  J»ZNW#A‘IW*T/ZW»TRRWtTRZW»TTTW»TZZP,TRRO,PHl 

**** 


TRANSFER  CAT*  ( GEOMETRIC)  OF  A~l  POINTS  FROM  TCNPQR ARY  STORAGE 
TO  PERMANENT  STORAGE 


»♦** 


3St  IF  I K- 1 »  353,3S3,362 
352  IF  UP.GT.23)  Gt  TO  bC'\ 

IF  (I  IP.GT.13I  .AND. I  IP.t.T.2'') )  GO  TO  6002 
RIAN! «J»«SR( Jl 
ZI  AMI  *  J  )•*  SZI  J  ) 

SR< J)=FAW 
SZIJI-ZNfc 

IF  IIIPLEF. GT.13I.ANO.IIPLEF.it.?'')>  go  TO  66*3 
GO  TO  6009 

IF  I(IPLfF.GT.13).AND.(lPLCF.LT.2P))  GO  TO  60-H 
GO  TO  b'VS 

6C04  T«IA*U.J)*TSR(  Jl 
TZIKMI,J»*TS7< J) 

<sOOR  TSRI  Jl  *RAW 
TSZI J)*ZAW 
GO  TO  60C9 

600?  IF  (IIPLFF.GT.13I.ANO.IIPLEF.LT.20) »  GO  TO  VP7 
6f  "  *  RTGHP*R  IKPli J) 

ZT  E 3P*  Z I A  FI . J  I 
R I  AMI »  J) *SRI J ) 

Z I  AMI «  J )*SZ I J) 

(33  SRI  J)  *PM» 

SZ( J) *ZNF 

oil  to 

6003  TRIAM1  ,  JJaTSRI J) 

TZIAM1,J)*TSZ<J) 
hCR  CONTINOF 

IF  (  IP.GT.23)  OF  Tfl  5P’6 
RHIKCOP, J)*«NF 
ZrtI  ACUP  t  J)*Z‘xH 
Ri)h(ACt)P,J)*«ONH 
ZOHIACUPf J)*ZPNH 
Go  TO  50 27 
5V6  TRH I A  f  J  )  *RNH 
TZ.OI  A ♦  J  )  *  ZNF 
TR  OH  I  A  »  J  )  *  Pn*jn 
TZOHIA.J)  *  Z^NF 
SV»  CciNT l MJF 

IF  I  l  TC-1  )  30  1 .1,  ’56,  3''0« 

3S  A  AIAC03*J)*Vw 
VI  AC  OR  t J) *VW 


P( ACORf  J  l=PW 
SZZIKCLF  .  J)  *S>Zv> 

SPRIKCUP  *  J)3SRR'» 
SRZI<CUR,J)=SRZA 
ST  TI  ACOR  *  J)  *STTt* 
TZZ<KCO»,J)='ZZW 
TPR ( ACUP  » J) *THPF 

TRZIACUR# J)*T«ZW 
TTTIACJP.J)=TTTW 
IF  IIP.G1.23)  cr  TO  300 
IF  ((  IP.CTol  3 1  .  ANO .  (  IP.LT.?'"' 


n0T  reproducible 

GO  TO  SO 30 
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woe 


r.o  r n  ir <r 
Hu  ru  *>031 

V'*  CiHTlMUF 

IP  ( IP.EC.l) 

1C  ALL  f  RAC  IKHl ,SHS,IPT  ,KUX  ,0  ,  Z ,  «H,  Z  H,  PITH,  ZW,  T»H  ,  T7  7  ,  Tk7  ,T  TT  , 

7  SPK,S77,SB/*SrT,P,  l^lll! 

IP  (1P.EC.U  CUl  FK AC TN  ITOR,T7Z,TBi,rTT,K,Kiix,P,?,Pnf7r»,kr>n,  »«»»«, 
IHUlf  Vim  tSPP  »S7».S«?,STT  ,P,IPT  } 

3'"*  CUNT1NUF 
KstMAX 

THIS  IS  THE  TO  CF  THE  K-lLUP 

**  **# 


f  i.tiW  K  IS  FCIAI  T (.  K«IU 

r 

^  ***** 


'c  j,<2-;nv,-7n 
tcin’-  ?NH-7u>» 

“1  O'* 

71  K(,  l«,J>  =SM  II 
RIkC'J«»J>  =  S*<  I) 

O.J  TO  I  er  3?  ,  37/*|  ,  IP 
•j'J*  ir  (  1PT(KMX,  I  l-l  *)  »?Af  37A.SA?3 

*»  )  A  *  T/(<MAX,J)  r  1',/U) 

TJ|K3AX,JJ  a  TS'(J) 

O'.  C.iMTlMlE 

if  (M.'c.n 

1C  ALL  FRAC(kmax,SHS.II»T»KiJX«K  ,  >  ,RH,7m.«I'm,  »nn, TKP , TZ Z , TP  7 ,  rrr, 

1  SKPf  S7Z,S«7,STT,p,  |V(ij, 

MJ  M  (  ’P",  OS  ,!•) 

>T,  1C  (  K  -  <  <4  <  ;j  )  »  )*,  iq",  VO 

)  7  •<  M.<txX(?! 

J.I  l-J'*!  7  > 

CALL  MPCVI  2,  srrxj7,TCCN7,  Kmm,»kk,  )JJ  I 

H  ’  r,i(xT|.vir 

!  '  Fuw  i AT  (thl,?iOn  ‘■«U,3y,MtTl^F=,^.At3H  *,»«.« 

]  VI  (A  ,3Xf7.<<»mi*,  IA//1  IX»P'|.'ELT  *,fP.4,xm  4S,’*, 

2  SlIL'P  I-  I HN  =  »*■".  i,3H  *S  ,'.K,thJ,K  ,  )  =  ,17,04) 

101  FORMAT  1  I  Sr  INTrUFAC*  *C»IL  *  ,  I  *  » 

1  2  FuP-lAT  I  7H  PF  C.l  v-N  »1  2  » 3*  »oHk  CCt  .  ,  l S// 3*  » IHJ  ,  11  x  ,  I  ■(/ ,  U *  ,  1  nu ,  S <  , 
:  TMiniz/  )*  i» < »  thsiiwfp  )  .sx,7  isii.i»7i  tSx»’<oi'..i  n  \ * , 

>  AH  SIC.  Ml  ,SX,',HSlr,<  O,0*»  3HPHJ/4H  ,  •'), 

.1  <■  (  SX,  ’Mix  U  ,  -5X  ,  ’r'Ov,  I 

M«PAT  I  !',,2F  1/  .4,  ,C1  2 .  ) ) 

3F  TJP  >. 

TN') 

SJttPUMTlLF  FpTFS')  UP*'n,?n,PA,7A,P3,/3,PC.;r,*w-  ,  .  ■  ’  . 

I ,  S R « w  .  S 7 / U  , S «  '  %  » S T T  w ,  x  , r .% p w ) 

1.4P11CIT  (.-CA|*q  (A-H,C-/,f),  INTT.,FP  (|-M 
N>M*n 


!>1  Z-\.  lMSO?-) 

‘jO  f't  lirlflfltlflf7|Pfip 

1 1  .  I  ,  l  .  1  ,  1  . 1  >  .  I  !» 

T/Zw=0.p 

TV//.: 

i>  V  =  -T  f  Tl*/  7 .0 
ST  T  4= -7  .r  *Pw 
SBZ<m.r 


io,i,i2,  i ,  n,  io,i,  i,  i , 


, 


T  .<1 
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m 


m 


W&r 


% 

$ 

i 

Ss 

I 


if 


I 


I 


t 


SZZW»PW 
GU  TU  1 
f  Z1«Z0 
RI*RO 
Z2«ZC 
R2«RC 
GO  TO  3C0 

1"  IF  tU.EC.2I.ANn.UR.F0.in  GO  TO  3 

zi*za 

Rl«RG 
R2*RA 
Z2*ZA 
R2*«A 
GO  TO  300 
12  Z1*ZA 
R1*RA 
Z2  *ZH 
R2*RR 
Gd  TO  TOO 
Z1  =ZB 
ai-RB 
Z2*ZC 
R?*OC 
GO  TO  200 

300  OUR1-02 

IF  tOABStOU-l. no-05)  301,361,302 

301  ETA*0.0 
GU  TU  310 

30^  D?*Z1-Z? 

IF  (OABS(rZ)-l. 00-05)  303,303,304 

’>•3  £1  A=P  l '  /?  .0 
GO  TO  3in 

3->4  ETA*0ATAM(Rl-R7)/<  Z1-Z2)  ) 

310  ALPHA=2.0*FTA 

TST*(a(  (  TRRV**T7ZW)  /2.  0  )■*•''.  5*  (  TZ  7W-T RRW )  *T)CUS  (  ALPHA) *TR Zw  • 
lllSINt  AlPFA) 

T«RW*n.5*TSTW*< 1 .r-0C0S( ALPHA) ) 

TZZW*->.5*TSTV»*t  l  .'HpCCSt  ALPHA)  ) 

TRZa=0.5«TSTW*OS IN) ALPHA) 

PW=-<  TRRM-T/ZV.UTTNI/3." 

SRRW*  TRRL+PV, 

SZZa«TZ3W*PW 
SRZW=TRZV 
STT\J»TTU*PH 
!  CONTINUE 


RETURN 

ENO 

/  /GO.  SYS  I N  on,  * 

1 

6031 

2  l 

10C0 
ICE 

i'.-' 

I 

AIM.) 

75000.7 

l  O.n 

/* 


NOT 


repp 


oouc'O^ 
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31 


r.?g  P73.onfl.rt 

0.7  C.o 

0.2R  1000000.0 


-7.7 


0.5 


0.Q2 

45 

2.7 


140.0 

?? 

UOOOO.O 

4b 


n.l?5 


n.125 


'.125 


n.125 
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